THERMODYNAMIC  MODELING  OF  PH  IN  CITRUS  SIMULATION  SYSTEMS 
UNDER  HIGH  PRESSURE  CARBON  DIOXIDE 


By 


BEHROUZ  MEYSSAMI 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 
OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 

UNIVERSITY  OF  FLORIDA 

1992 


UNIVERSITY  OF  FLORIDA  IffO 


To  my  wife,  parents  and  brother 


i 


ACKNOWLEDGEMENTS 

I  would  like  to  express  my  sincere  gratitude  to  Dr. 
Arthur  A.  Teixeira  and  Dr.  Murat  0.  Balaban,  my  major 
advisors,  for  their  support,  guidance  and  assistance 
throughout  the  course  of  this  study.  I  am  also  very  grateful 
to  Dr.  Khe  V.  Chau,  Dr.  Robert  J.  Braddock  and  Dr.  Spyros  A. 
Svoronos,  members  of  my  advisory  committee,  for  their 
guidance,  suggestions  and  assistance. 

I  am  thankful  to  Linda  L.  Latta,  our  laboratory 
technician  for  her  suggestions  and  assistance  in  working  with 
the  equipment. 

I  am  very  grateful  to  my  wife,  Kimia,  my  parents, 
Mohtaram  and  Asdollah  Meyssami  and  my  brother,  Bahram,  for  all 
the  support  and  encouragement  they  have  given  me  in  my 
studies. 


iii 


TABLE  OF  CONTENTS 


ACKNOWLEDGEMENTS    iii 

LIST  OF  TABLES   vii 

LIST  OF  FIGURES   ix 

KEY  TO  SYMBOLS  xi 

ABSTRACT   xiv 

CHAPTERS 

1  INTRODUCTION    1 

2  REVIEW  OF  LITERATURE    4 

Supercritical  C02    4 

Properties  and  Applications    4 

Thermodynamic  Models  of  Phase  Equilibria     ...  7 

Equations  of  State    9 

Liquid  Supercritical  Phase  Equilibrium  Models  10 
Pectinase  Action  and  Citrus  Juice  Cloud  ....  12 
Pectic  Substances  and  Pectic  Enzymes     ....  12 

Classification  of  Pectic  Enzymes    13 

Mechanism  of  Action  of  Pectinases    14 

Composition  of  Citrus  Juices    15 

Modeling  of  pH  Using  Reaction  Equilibria  in  the 

Simulation  Model  systems    16 

Electrolyte  Solutions  and  Their  Modeling  .  .  16 
Equilibrium  Reactions  in  the  Liquid  Phase  .  .  18 
Effect  of  Dissolved  Carbon  Dioxide  on  pH  of 

Pure  Water   18 

Prediction  of  Thermodynamic  Properties  of 
Aqueous  Electrolytes  at  High  Pressures     ...  26 
Pressure  Dependence  of  Molar  Volume  ....  26 
Pressure  Dependence  of  Equilibrium 

Constant   28 

Pressure  Dependence  of  Activity 

Coefficient  3  0 


iv 


3     MATERIALS  AND  METHODS 


32 


Development  of  Thermodynamic  Models    32 

Calculation  of  Activity  Coefficients    32 

Pierotti  et  al.  Method  for  Carbonic  Acid     .   .  32 

Debye-Huckel  Method  for  Individual  Ions  ...  34 
UNIFAC  Method  for  Ascorbic  and  Citric 

Acids   36 

Calculation  of  pH  of  the  Simulation  Systems  .  .  42 
Systematic  Method  for  Calculating  Equilibrium 

Concentrations    43 

Ionization  Constant  of  Water    43 

Numerical  Values  of  Equilibrium  Constants  .  .44 
Set  of  Equations  for  the  Pure  Water-Carbon 

Dioxide  System    45 

Set  of  Equations  for  the  Ascorbic  Acid- 
Water-Carbon  Dioxide  System    48 

Set  of  Equations  for  the  Citric  Acid- 
Water-Carbon  Dioxide  System    49 

Set  of  Equations  for  the  Ascorbic  Acid- 
Citric  Acid-Water-Carbon  Dioxide  System  .  .51 

Experimental    53 

High  Pressure  pH  Probe   53 

Experimental  Design    54 

Preliminary  Experiments    55 

Experiments  on  the  Effect  of  Pressure  and  Acid 

Concentrations  on  the  pH  of  Model  Solutions  .  55 

Statistical  and  Sensitivity  Analyses    60 

4     RESULTS  AND  DISCUSSION    62 

Pure  Water-Carbon  Dioxide  Simulation  System  ...  62 

Calculated  Activity  Coefficients    62 

Variation  of  Activity  Coefficients  with 

Pressure   63 

Variation  of  Equilibrium  Constants  with 

Pressure   67 

Comparison  of  the  Experimental  and  Calculated  pH 

Values   68 

Effect  of  Extreme  Process  Pressures    68 

Effect  of  Sugar  on  Solubility  of  Carbon 

Dioxide   71 

Effect  of  process  temperature  on  pH   74 

Ascorbic  Acid-Water-Carbon  Dioxide  Simulation 

System   74 

Calculated  Activity  Coefficients    74 

Comparison  of  Calculated  and  Experimental  pH 

Values   77 

Citric  Acid-Water-Carbon  Dioxide  Simulation 

System   82 

Calculated  Activity  Coefficients    82 


v 


Comparison  of  Calculated  and  Experimental  pH 

Values  83 

Ascorbic  Acid-Citric  Acid-Water-Carbon  Dioxide 

Simulation  System    88 

Calculated  Activity  Coefficients    88 

Comparison  of  Calculated  and  Experimental  pH 

Values  88 

Single  Strength  Orange  Juice    97 

Model  Sensitivity  Analysis    97 

Discussion  of  Results    101 

5     CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  STUDY     .  105 

Conclusions  105 

Suggestions  For  Future  Study    107 

APPENDIX  A       TABLES  OF  PREDICTED  AND  MEASURED  PH  VALUES  108 

APPENDIX  B       FORTRAN  PROGRAM  TO  SOLVE  FOR  EQUILIBRIUM 

ACTIVITIES  USING  NEWTON-RAPHSON  ITERATION 
METHOD  115 

REFERENCE  LIST  126 

BIOGRAPHICAL  SKETCH  13  6 


vi 


LIST  OF  TABLES 


Table  Page 

3-1  Pierotti  parameters  for  calculating  the  infinite 
dilution  activity  coefficient  of  carbonic  acid 
in  water  34 

3-2    Approximate  effective  ionic  radii  of  inorganic  and 

organic  ions  35 

3-3    Group  volume  and  surface-area  parameters  for 

citric  and  ascorbic  acids  used  in  the  UNIFAC 

method  39 

3-4    Ionization  constant  for  water  (KJ  44 

3-5    Dissociation  constants  for  carbonic,  ascorbic  and 

citric  acids    45 

3-  6    Coded  design  variable  settings  for  the  pH 

experiments  56 

4-  1    Calculated  activity  coefficients  for  water-carbon 

dioxide  system    65 

4-2    Compressibility  coefficients  of  the  ions  in  the 

water-carbon  dioxide  system    65 

4-3    Measured  pH  values  for  pure  water  and  sucrose 

solution  75 

4-4    Measured  pH  values  of  water-carbon  dioxide  system 

at  different  process  temperatures    76 

4-5    Calculated  activity  coefficients  for  the  ascorbic 

acid-water-carbon  dioxide  system    77 

4-6    Calculated  activity  coefficients  for  the  citric 

acid-water-carbon  dioxide  system    83 

4-7    Sensitivity  analysis  of  the  pH  model  at 

atmospheric  pressure  and  37°C  with  the  variation 

in  activity  coefficients    103 

vii 


4-8    Sensitivity  analysis  of  the  pH  model  at 

atmospheric  pressure  and  37°C  with  the  variation 

in  the  equilibrium  constants  104 

A-l  Predicted  and  measured  pH  values  at  different 
process  pressures  for  the  pure  water-carbon 
dioxide  simulation  system    109 

A-2    Comparison  of  measured  and  predicted  pH  values  of 
the  ascorbic  acid-water-carbon  dioxide  system 
containing  1.107x10"   M  ascorbic  acid    no 

A-3    Comparison  of  measured  and  predicted  pH  values  of 
the  ascorbic  acid-water-carbon  dioxide  system 
containing  5.54x10"   M  ascorbic  acid    110 

A-4    Comparison  of  measured  and  predicted  pH  for  the 
citric  acid-water-carbon  dioxide  system  with 
6.08x10"   M  citric  acid   Ill 

A-5    Comparison  of  measured  and  predicted  pH  for  the 
citric  <icid-water-carbon  dioxide  system  with 
3.04x10"   M  citric  acid   Ill 

A-6    Comparison  of  measured  and  predicted  pH  for  the 
ascorbic  acid-citric  acid-water-carbon  dioxide 
system  with  5.54x10"    M  ascorbic  acid  and 
6.08x10"   M  citric  acid   112 

A-7    Comparison  of  measured  and  predicted  pH  for  the 
ascorbic  acid-citric  acid-water-carbon  dioxide 
system  v^ith  1.107x10"   M  ascorbic  acid  and 
6.08x10    M  citric  acid   112 

A-8    Comparison  of  measured  and  predicted  pH  for 

ascorbic  acid-citric  acid-water-carbon  dioxide 

system  with  1.107x10"    M  ascorbic  acid  and 

3.04x10    M  citric  acid   113 

A-9    Comparison  of  measured  and  predicted  pH  for 

ascorbic  acid-citric  acid-water-carbon  dioxide 

system  with  5.54x10"   M  ascorbic  acid  and 

3.04x10    M  citric  acid   113 

A-10  Comparison  of  measured  and  predicted  pH  of  single 

strength  orange  juice    114 


viii 


LIST  OF  FIGURES 

Figure  Page 

2-1    Liquid-phase  chemical  reactions  in  a  water-carbon 

dioxide  binary  system    20 

2-2    Dissociation  reactions  of  ascorbic  acid    22 

2-  3    Dissociation  reactions  of  citric  acid    24 

3-  1    Schematic  of  the  experimental  setup  for  measuring 

pH  in  the  supercritical  vessel  at  high 

pressures  59 


4-1    Measured  and  predicted  pH  of  water-carbon  dioxide 

simulation  system  at  pressures  up  to  34.48  MPa  .   .  70 

4-2    Measured  and  predicted  pH  of  water-carbon  dioxide 

simulation  system  at  pressures  upto  3.363  MPa     .   .  73 

4-3    Measured  and  predicted  pH  of  the  ascorbic  acid- 


water-car^on  dioxide  system  containing 
1.107x10"  M  ascorbic  acid  (HI:A-L0W:C)    79 

4-4    Measured  and  predicted  pH  of  the  ascorbic  acid- 
water-car^bon  dioxide  system  containing 
5.540x10"  M  ascorbic  acid  (MID:A-L0W:C)   81 


4-5    Measured  and  predicted  pH  of  the  citric  acid- 
water-ca^rbon  dioxide  system  containing 
6.08x10"   M  citric  acid  (LOW:A-HI:C)  .... 

4-6    Measured  and  predicted  pH  of  the  citric  acid- 
water-ca^rbon  dioxide  system  containing 
3.04x10    M  citric  acid  (L0W:A-MID:C)       .   .  . 

4-7    Measured  and  predicted  pH  of  the  ascorbic  acid- 
citric  acid-water-carbon  dioxide  system 
containing  5.54x10"   M  ascorbic  acid  and 
6.08x10    M  citric  acid  (MID:A-HI:C)  .... 

4-8    Measured  and  predicted  pH  of  the  ascorbic  acid- 
citric  acid-water-carbon  dioxide  system 
containing  1.107x10"   M  ascorbic  acid  and 
6.08x10     M  citric  acid  (HI:A-HI:C)  .... 


ix 


4-9    Measured  and  predicted  pH  of  the  ascorbic  acid- 
citric  acid-water-carbon  dioxide  system 
containing  1.107x10*   M  ascorbic  acid  and 
3.04x10    M  citric  acid  (HI:A-MID:C)    94 

4-10  Measured  and  predicted  pH  of  the  ascorbic  acid- 
citric  acid-water-carbon  dioxide  system 
containing  5.54x10"   M  ascorbic  acid  and 
3.04x10    M  citric  acid  (MID: A-MID: C)    96 

4-11  Measured  and  predicted  pH  of  single  strength 

orange  juice    99 


x 


KEY  TO  SYMBOLS 


SYMBOL  DEFINITION  UNITS 

M  Chemical  potential  J/mole 

A  Helmholtz  free  energy  J 

n  Moles  moles 

T  Temperature  °C 

Y  Molar  volume  cm3/mole 
V°  Standard  molar  volume  cm3/mole 
P                       Process  pressure  MPa 

P°  Standard  pressure  MPa 

AV  Change  of  molar  volume  cm3/mole 

AVR  Change  of  reaction  volume  cm3 

k°  Standard  isothermal 

compressibility  coefficient  cm3/mole.KPa 

v  Stoichiometric  coefficient  (dimensionless) 

K  Eguilibrium  constant  moles/cm3 

a  Activity  moles/cm3 

Y  Activity  coefficient  (dimensionless) 

c 

Y  Combinatorial  activity 

coefficient  (dimensionless) 

Y  Residual  activity 

coefficient  (dimensionless) 


xi 


Molar  concentration 

Ideal  gas  constant 

Ionic  radius 

Ionic  strength 

Ionic  charge 

Mole  fraction 

Group  numbers 

van  der  Waals  volume 
parameter 

van  der  Waals  surface 
parameter 

Group  surface  area 
fraction 

Group  fraction 

Group  interaction 
parameter 

Exponential  of  group 
interaction  parameter 
over  temperature 

Solubility 

Coded  design  variable 

Design  variable 

Mean  of  design  variables 

Scale  factor 

Test  statistic 

Null  hypothesis 

Alternative  hypothesis 

Average  difference  between 
treatment  means 


moles/cm 
J/g.mole.K 
Angstroms 
moles/cm3 
(dimensionless) 
(dimensionless) 
(dimensionless) 

(dimensionless) 

(dimensionless) 

(dimensionless) 
(dimensionless) 

J 

exp(J/K) 
moles/cm3 
(dimensionless) 


xii 


Difference  between  treatment 

means   

Standard  deviation  of  paired 

differences   

Coefficient  of  correlation  (dimensionless) 


xiii 


Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

THERMODYNAMIC  MODELING  OF  PH  IN  CITRUS  SIMULATION  SYSTEMS 
UNDER  HIGH  PRESSURE  CARBON  DIOXIDE 

By 

Behrouz  Meyssami 
May  1992 

Chairman:  A.  A.  Teixeira 
Cochairman  :  M.  O.  Balaban 

Major  Department:  Agricultural  Engineering 

Thermodynamic  models  were  developed  to  predict  the  pH 
of  liquid  citrus  simulation  systems  as  functions  of  high 
pressure  carbon  dioxide  process  conditions.  Such  processing 
holds  promise  as  an  alternative  to  heat  pasteurization  for 
inactivating  pectic  enzymes  through  temporarily  lowering  pH 
instead  of  using  heat.     Predicted  pH  values  were  compared 
with  those  measured  experimentally  at  different  process 
conditions  with  the  aid  of  a  pH  probe    especially  designed 
to  withstand  high  process  pressures.  The  four  simulation 
systems  included  pure  water-carbon  dioxide,  ascorbic  acid- 
water-carbon  dioxide,  citric  acid-water-carbon  dioxide  and 
ascorbic  acid-citric  acid-water-carbon  dioxide.  Gage 
process  pressures  ranged  from  0  to  34.50  MPa  at  temperatures 
of  32,  37  and  42°C.     Acid  concentrations  ranged  from  0  to 
1.107xl0"3  M  and  0  to  6.080xl0"4  M  for  ascorbic  and  citric 
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acids,  respectively.     For  all  the  simulation  systems  the 
predicted  pH  values  closely  followed  the  measured  pH  values. 
In  the  pure  water-carbon  dioxide  system,  the  pH  was  a  strong 
function  of  process  pressure  but  a  weak  function  of  process 
temperature.     In  the  ascorbic  acid-water-carbon  dioxide, 
citric  acid-water-carbon  dioxide  and  ascorbic  acid-citric 
acid-water-carbon  dioxide  systems,  however,  the  pH  was 
relatively  insensitive  to  pressure  while  remaining  a  strong 
function  of  acid  concentrations. 

The  pH  in  the  ternary  and  quaternary  simulation  systems 
showed  a  very  small  variation  over  the  experimental  pressure 
range,  indicating  that  low  pH  alone  can  not  totally  account 
for  pectinesterase  inactivation  in  citrus  juice  treated  with 
supercritical  carbon  dioxide.     Since  the  pH  variation  was 
greater  in  the  water-carbon  dioxide  system  than  the  other 
three  systems,  more  carbon  dioxide  was  dissolved  in  the 
binary  system  resulting  in  the  generation  of  greater 
quantities  of  carbonic  acid  in  the  solution.     The  reduction 
in  the  solubility  of  carbon  dioxide  in  the  solution  was 
believed  to  be  caused  by  the  presence  of  other  ions  that 
formed  as  a  result  of  dissociation  of  ascorbic  and  citric 
acids  in  the  ternary  and  guaternary  systems. 
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CHAPTER  1 
INTRODUCTION 


The  cloud  in  citrus  juice  is  responsible  for  most  of 
the  characteristic  color  and  appearance  that  is  widely 
recognized  and  appreciated  by  consumers.     The  retention  of 
this  cloud  in  most  citrus  juices,  concentrates  and  beverages 
is  of  prime  importance  to  the  citrus-processing  industry  for 
consumer  acceptance.  Most  of  the  published  work  relating  to 
the  cloud  loss  in  citrus  juice  products  attributes  the  loss 
mainly  to  the  activity  of  the  pectic  enzymes,  particularly 
pectinesterase  naturally  present  in  the  juice  (Baker,  1977) . 
Part  of  the  methylester  linkages  of  the  soluble  pectin  is 
split  by  pectinesterase,  resulting  in  the  formation  of  low- 
methoxyl  pectin  which,  in  the  presence  of  calcium  and 
magnesium  cations,  precipitates  and  results  in  the  cloud 
loss.  Divalent  cations  must  be  present  in  solution  for  the 
pectin  molecules  to  agglomerate  and  subsequently  settle  out 
(Rouse  and  Atkins,  1952) . 

Most  procedures  for  prevention  of  cloud  loss  are  based 
on  heat  inactivation  of  pectinesterase  by  pasteurization 
(Rouse  and  Atkins,  1952;  Rouse  and  Atkins,  1953;  Mizrahi  and 
Berk,  1970;  Versteeg  et  al.,  1980;  Wicker  and  Temelli, 
1988) .  However,  pasteurization    temperatures  can  result  in 
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flavor  and  aroma  changes,  as  well  as  color  changes  due  to 
nonenzymatic  browning. 

Owusu-Yaw  et  al.   (1988)  have  shown  that  pectinesterase 
can  also  be  inactivated  by  lowering  the  pH  of  the  juice 
using  either  a  cation  exchange  resin  or  by  direct  addition 
of  hydrochloric  acid  to  the  orange  juice.  Treatment  with 
supercritical  carbon  dioxide  has  been  recently  proposed  as  a 
promising  alternative  process  for  inactivating 
pectinesterase  and  improving  cloud  stability  (Arreola, 
1990) .  Cloud  stability  in  this  process  is  based  on  the 
hypothesis  that  under  pressure  carbon  dioxide  dissolves  in 
the  solution  producing  carbonic  acid,  thereby  temporarily 
lowering  the  pH  and  inactivating  pectinesterase.  Arreola 
(1990)   showed  that  pectinesterase  was  inactivated  in  the 
supercritical  process  with  various  combinations  of  pressure, 
temperature  and  process  time. 

Studying  the  variation  of  pH  during  the  supercritical 
process  provides  the  means  of  better  understanding  how 
pectinesterase  is  inactivated  in  this  process.  Also,  for 
carbonated  beverages  it  is  important  to  know  what  the  pH 
variation  is  during  the  carbonation  process  as  well  as  when 
it  is  packed  under  pressure  (Moore  and  Buchanan,   1930) . 

Citrus  juices  and  beverages,  with  their  complex 
composition,  can  be  viewed  as  electrolyte  solutions  of 
multiple  ions.  During  the  supercritical  process,  the  pH 
varies  with  the  process  pressure,  temperature  and  solution 
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composition.  From  the  chemical  process  design  point  of  view, 
the  studies  of  food  electrolyte  solutions  and  applications 
of  the  supercritical  carbon  dioxide  process  to  these  systems 
have  rarely  been  done  in  the  past  (Maurer,   1983) .  In  the 
supercritical  carbon  dioxide  process,  through  the  use  of 
reaction  equilibria,  the  pH  of  the  solution  may  actually  be 
calculated  at  specific  process  conditions.  This  involves 
finding  a  mathematical  expression  (model)  for  the  pH  as  a 
function  of  the  process  conditions  using  thermodynamic 
relations  and  available  data. 

The  objectives  of  this  study  were  to 

(1)  Develop  thermodynamic  models  to  predict  the  pH  of 
four  citrus  simulation  systems  of  (a)  pure  water-carbon 
dioxide  ,     (b)  ascorbic  acid-water-carbon  dioxide, 

(c)  citric  acid-water-carbon  dioxide  and  (d)  ascorbic  acid- 
citric  acid-water-carbon  dioxide. 

(2)  Test  the  models  by  comparing  the  predicted  pH  with 
the  actual  measured  pH  values  for  each  system  at  different 
process  conditions. 

(3)  Determine  significance  of  pH  lowering  in 
pectinesterase  inactivation  and  subsequent  citrus  cloud 
stability  by  studying  variations  in  pH  over  the  experimental 
pressure  range  at  different  temperatures  and  system 
compositions. 


CHAPTER  2 
REVIEW  OF  LITERATURE 

2 . 1  Supercritical  CO, 

2.1.1  Properties  and  Applications 

Supercritical  carbon  dioxide  processing  is  a  unique 
unit  operation  that  utilizes  the  special  properties  of  the 
solvent  carbon  dioxide  above  its  critical  temperature  and 
pressure.  The  development  of  the  supercritical  process 
technology  has  largely  been  based  on  practical  experience 
and  engineering  intuition.  The  process  can  perform  different 
operations,  depending  on  the  nature  of  the  starting 
material,  the  operating  conditions,  and  the  supercritical 
solvent  (Rizvi  et  al.,  1986a).  Each  of  these  factors  must  be 
carefully  considered  before  proceeding  with  the  development 
of  a  successful  supercritical  process. 

Carbon  dioxide  has  a  critical  temperature  of  31.1°C  and 
a  critical  pressure  of  7.38  MPa  (McHugh  and  Krukonis,  1986). 
Within  the  supercritical  region  carbon  dioxide  exhibits  its 
special  properties  such  as  low  viscosity,  low  surface 
tension,  and  high  diffusivity  (de  Fillippi,  1982).  In  this 
region,  carbon  dioxide  exhibits  physicochemical  properties 
intermediate  between  those  of  liquids  and  gases,  which 
enhances  its  role  as  a  solvent  .  Its  relatively  high  density 
results  in  its  good  solvent  power,  while  its  relatively  low 
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viscosity  and  diffusivity    provide  appreciable  penetrating 
power  into  the  solute  matrix.  Chrastil  (1982)  experimentally 
observed  that  the  concentration  of  a  solute  dissolved  in  a 
supercritical  fluid  is  related  to  the  density  of  the  fluid, 
raised  to  some  positive  power.  This  observation,  coupled 
with  the  pressure-density-temperature  behavior  indicates 
that  the  greatest  solubility-pressure  and  solubility- 
temperature  variations  occur  in  the  region  near  the  critical 
point  of  the  solvent.  As  the  temperature  increases  within 
this  region,  the  density  of  the  solvent  decreases  rapidly, 
and  the  solubility  of  a  solute  would  be  expected  to  decline 
despite  its  increased  vapor  pressure.  At  very  high 
pressures,  however,  the  vapor  pressure  of  the  solute  varies 
more  strongly  with  temperature  than  the  density  of  the 
solvent.  Thus,  the  solubility  of  a  solute  increases  with 
increasing  temperature  at  high  pressure.  Similarly,  a 
decrease  in  the  solvent  power  of  a  compressed  gas  can  be 
expected  when  pressure  is  reduced.  These  rough  qualitative 
conclusions  from  the  density-pressure  behavior  of  a  solvent 
were  discussed  by  Kurnick  et  al.   (1980;  1981)  working  on  the 
solubility  of  benzoic  acid  in  carbon  dioxide. 

The  chemical  complexity  of  most  food  ingredients  and 
their  vulnerability  to  react  and  degrade  at  elevated 
temperatures  emphasize  the  considerations  that  must  be  taken 
into  account  in  supercritical  solvent  selection.  Because  of 
its  unique  properties,  carbon  dioxide  has  a  wide  range  of 
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food  processing  applications.  It  is  an  inert,  nonflammable, 
nontoxic  and  nonexplosive  substance  with  a  moderate  critical 
pressure  which  minimizes  compression  costs.  Among  the  many 
applications  of  supercritical  carbon  dioxide,  extraction  and 
fractionation  cover  a  good  percentage.  They  include 
decaf feination  of  coffee,  recovery  of  hop  extracts  and  oil 
and  fat  removal  applications  (Bott,  1982;  Vollbrecht,  1982; 
Calame  and  Steiner,  1982;  Gardner,  1982). 

Beyond  extraction  and  fractionation,  carbon  dioxide  has 
considerable  antimicrobial  activity.  The  degree  of 
inhibition,  as  reported  by  Clark  and  Takacs  (1980) ,  depends 
on  the  type  of  microorganism.  Pseudomonas  appear  to  be  very 
sensitive  while  other  types  of  organisms  such  as 
lactobacilli,  yeasts,  and  Clostridia  are  less  sensitive 
(Jones  and  Greenfield,   1982;  Molin,  1983).  Carbon  dioxide 
also  acts  synergistically  with  certain  antimicrobial  agents 
such  as  sorbate  (Haas  and  Herman,  1977a)  and  propylene 
glycol  (Haas  and  Herman,  1977b) . 

The  inhibitory  effect  of  carbon  dioxide  has  also  been 
reported  to  increase  when  it  is  applied  under  pressure 
(Kauffman  et  al.,  1969;  Enfors  and  Molin,  1978,  1980; 
Blickstad  et  al.,   1981;  Doyle,   1983;  Kamihira  et  al.,  1987). 
Haas  et  al.   (1989)  determined  that  the  moisture  content  (or 
water  activity,  aH)  of  the  product  was  very  important  to  the 
antimicrobial  action  of  C02  in  wheat  flour,  while  the 
temperature  of  exposure  was  not.  There  was  slightly  less 
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effect  at  -22°C  than  at  23°C  and  50°C.  Haas  et  al.  (1989) 
also  performed  experiments  with  freshly  squeezed  orange 
juice  (pineapple  variety)  and  reported  that  the  microflora 
consisting  of  gluconobacter ,  yeast,  mold,  and  lactic  acid 
bacteria  were  completely  eliminated  at  a  pressure  of  5.516 
KPa.  Although  not  much  data  are  available  on  the 
inactivation  of  pectinases  by  carbon  dioxide,  Haas  et  al. 
(1989)  reported  a  24-48  hour  delay  in  flocculation  of  orange 
juice  cloud.  They  speculated  that  either  pectinases  may  be 
inactivated  at  a  pressure  of  5.516  KPa  or  it  could  be  that 
some  of  the  C02  remains  bound  to  the  cloud  affecting  its 
specific  gravity. 

2.1.2  Thermodynamic  Models  of  Phase  Equilibria 

Mixtures  of  materials  at  supercritical  conditions 
exhibit  highly  nonideal  behavior  and  do  not  lend  themselves 
easily  to  quantitative  data  correlation  and  phase 
equilibrium  predictions.  This  is  particularly  true  for 
materials  of  biological  origin,  because  their  physical 
property  data  are  usually  extremely  scarce —  if  available  at 
all  (Rizvi  et  al.,  1986b). 

Models  of  equilibrium  between  supercritical  and 
condensed  phases  can  be  separated  into  two  general  classes: 
those  describing  supercritical  fluid-liquid  equilibria  and 
those  describing  supercritical  fluid-solid  equilibria.  As 
with  all  other  equilibria,  the  condition  for  equilibrium  is 
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that  the  chemical  potential  (/i,)  of  each  component  i  in  one 
phase  be  equal  to  the  chemical  potential  of  that  component 
in  all  other  phases,  that  is, 

11/  -  |li>  -  (2-D 
where  i=l,2,...,N  for  an  N-component  system  in  phases 

(2,3,0),  .  •  • 

The  main  differences  between  models  of  supercritical 
fluid-liquid  and  supercritical  fluid-solid  mixtures  are  that 
the  solid  phase  is  assumed  to  dissolve  no  solvent 
(Prausnitz,  1969;  Gitterman  and  Procaccia,  1983;  Procaccia 
and  Gitterman,  1983) ,  and  in  multicomponent  systems,  no 
solid  solutions  are  assumed  to  form  (Kurnick  and  Reid, 
1982) .  These  two  assumptions  simplify  the  analysis  by 
eliminating  the  composition  dependence  of  the  chemical 
potential  of  the  solid  phase.  In  case  of  equilibrium  between 
liquids  and  supercritical  fluids,  however,  multicomponent 
solutions  are  common  and  considerable  amounts  of 
supercritical  solvent  can  dissolve  into  the  liquid  phase 
(Brunner  and  Peter,  1982)  leading  to  the  composition 
dependence  of  the  chemical  potential. 


2.1.3  Equations  of  State 

Phase  equilibrium  calculations  require  a  method  for 
computing  the  chemical  potential  from  available  data.  The 
chemical  potential  of  any  component  i  in  a  supercritical 
fluid  or  liquid  phase  a  may  be  computed  from  the  molar 
Helmholtz  free  energy  A  by  means  of  the  Dieters  and 
Schneider  (1976)  relation: 


where  the  Helmholtz  free  energy  is  a  function  only  of  the 
pressure,  volume,  temperature  and  composition.  Knowing  the 
Helmholtz  free  energy  at  the  reference  state  of  a  pure 
component  and  the  chemical  potential  (calculated  by  an 
equation  of  state) ,  the  Helmholtz  energy  can  then  be 
predicted  at  the  equilibrium  state.  Numerous  process- 
oriented  equations  of  state  which  describe  high  density 
fluids  have  been  proposed.  These  equations  vary  in  accuracy, 
number  of  parameters,  and  amount  and  type  of  data  required 
to  fix  their  parameters.  Vidal  (1983)  and  Gubbins  (1983) 
offer  good  reviews  of  equations  of  state  and  mixing  rules 
for  their  parameters.  Pure-component  equation  of  state 
parameters  for  materials  for  which  pressure-volume- 
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temperature  are  unavailable  may  be  back-calculated  from 
phase  equilibrium  data.  Mixture  parameters  are  always  back- 
calculated  from  binary  or  multicomponent  phase  equilibria 
data  usinq  nonlinear  reduction  techniques  (Anderson  et  al., 
1978,  and  Skjold-Jorqensen,   1983).  Some  of  the  commonly  used 
equations  of  state  include  van  der  Waals,  Redlich-Kwonq 

(1948) ,  Penq-Robinson  (1976) ,  and  Teja-Patel  equations 

(Rizvi  et  al.,  1986b). 

2.1.4  Liquid-Supercritical  Phase  Equilibrium  Models 

Models  of  liquid-supercritical  phase  equilibria  can  be 
classified  into  two  types:  those  which  treat  the  liquid 
phase  as  a  dense  qas,  and  those  which  treat  the  system  as  a 
classical  liquid-vapor  equilibrium.  The  first  approach  has 
been  shown  to  be  useful  by  van  Konyenburq  (1968) ,  who  proved 
that  the  van  der  Waals  equation  of  state  could  predict  five 
out  of  the  six  known  classes  of  fluid  phase  diaqrams  and 
predicted  phase  behavior  that  had  not  yet  been 
experimentally  observed.  In  this  approach,  an  equation  of 
state  is  used  to  describe  both  the  liquid  and  supercritical 
phases,  by  computinq  the  chemical  potential  from  an  equation 
of  state  for  each  component  in  the  system.  To  apply  these 
equations  to  biomaterials  for  which  equations  of  state  are 
unknown  and  critical  data  are  not  available,  binary 
solubility  data  are  used  to  obtain  estimates  of  the  equation 
of  state  parameters  for  the  nonvolatile  solute.  Solute- 
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solvent  interaction  parameter (s)  are  estimated  from 
solubility  data,  even  if  the  equation  of  state  parameters 
are  unknown  for  the  solute.  To  extend  the  calculation  to 
ternary  systems,  the  parameters  obtained  in  the  way  just 
described,  from  binary  data,  may  be  used  in  conjunction  with 
some  ternary  data  to  estimate  the  interaction  parameters 
(Rizvi  et  al.,  1986b). 

This  model  has  been  shown  to  adequately  represent 
vapor-liquid  equilibria  for  numerous  systems  of 
hydrocarbons  and  low  molecular  weight  molecules  over  a  wide 
range  in  polarity  and  sometimes  even  close  to  critical 
points  (Gray  et  al.,  1970;  Soave,  1972;  Patel  and  Teja, 
1982;  Mollerup  and  Fredenslund,  1979).  King  et  al.  (1983) 
have  used  a  modified  Soave-Redlich-Kwong  equation  of  state 
for  applying  the  model  to  various  binary  systems  involving 
triglycerides  mixtures  and  fatty  acids.  Their  results 
indicate  generally  good  correlation  between  predicted  and 
experimental  phase  equilibrium  data,  although  precise 
experimental  data  is  difficult  to  obtain  in  the 
supercritical  region. 

Correlations  based  on  vapor  pressure  data  (Joffe  et 
al.,  1970)  and  liquid  molar  volume     (Brunner  and  Hederer, 
1979)  have  been  developed.  These  correlations  permit 
relatively  precise  estimation  of  Redlich-Kwong  pure- 
component  equation  of  state  parameters  for  heavy 
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biomaterials,  such  as  long-chain  fatty  acids  and 
triglycerides . 

In  the  classical  vapor-liquid  equilibrium  model, 
activity  coefficients  are  used  to  correlate  experimental 
data.  Since  activity  coefficients  are  not  required  to  have 
any  compositional  functionality,  the  method  has  the 
advantage  that  data  which  may  not  be  used  with  an  equation 
of  state  can  be  used  with  equations  such  as  Henry's  law 
equation  to  predict  the  behavior  of  the  components 
(Prausnitz,   1969;  Reid  et  al . ,   1977).   From  this  point  of 
view,  if  enough  physical  property  data    as  well  as  activity 
coefficient  data  could  be  obtained  from  literature,  this 
technique  would  provide  a  powerful  and  fundamentally  valid 
way  to  correlate  data,  especially  if  the  data  cannot  be 
adequately  correlated  by  equations  of  state. 

2.2  Pectinase  Action  and  Citrus  Juice  Cloud 
2.2.1  Pectic  Substances  and  Pectic  Enzymes 

The  pectic  substances  in  the  plant  tissues  and  the 
enzymes  which  degrade  them  are  of  major  importance  to  the 
food  industry  because  of  their  effect  on  texture  of  food  and 
in  the  preparation  of  wines  and  fruit  juices.  The 
application  of  enzymes  in  the  production  of  fruit  juices, 
extracts  and  concentrates,  was  first  introduced  by  Kertesz 
(1930,1931),  Willaman  (1933)  and  Willaman  and  Kertesz  (1931) 
in  the  United  States.  Although  many  articles  have  been 


published  in  this  area,  the  mechanism  of  action  of 
pectinases  is  not  fully  understood.  Kertesz  (1930)  reported 
that  the  cloudiness  of  fruit  juices  was  partly  due  to  their 
pectin  content  where  the  colloidally  suspended  pectin  would 
hold  other  substances  in  suspension.  Pectic  enzymes 
(pectinases)  have  been  designated  by  Kertesz  as  that  group 
of  enzymes  which  hydrolyze  pectic  substancess  into  reducing 
sugars.  Pectic  (Greek=congealing)  substances  are  a  group  of 
complex  colloidal  plant  carbohydrates  made  up  primarily  of 
Alpha- (1-4) -D-polygalacturonic  acid  located  mostly  in  the 
middle  lamella  between  cells  in  higher  plant  tissues. 
Polygalacturonic  acid  exists  in  varying  states  of 
methoxylation  and  neutralization.  Pectin  is  composed  of 
polygalacturonic  acid  with  about  75%  of  the  carboxyl  groups 
being  esterified  with  methanol.  Pectic  acid  is  an  almost 
pure,  colloidal,  polygalacturonic  acid  free  from  methoxyl 
groups.  Normal  or  acid  salts  of  pectic  acid  are  called 
pectates  (Willaman,  1933) . 

2.2.2  Classification  of  Pectic  Enzvmes 

Pectinases  can  be  divided  into  two  categories: 
Depolymerizing  enzymes  and  saponifying  enzymes  (or 
pectinesterases,  P.E.).  The  depolymerizing  enzymes  are 
further  classified  according  to:  hydrolytic  or 
transeliminative  cleavage  of  glycosidic  bonds,  endo-  or  exo- 
mechanism  of  the  splitting  action,  i.e.,  whether  breakdown 
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starts  from  within  the  polymer  or  from  the  end  of  the 
polymer;  and  preference  for  pectic  acid  or  pectin  as  a 
substrate.  Both  polymethylgalacturonase  (PMG)  and 
polygalacturonase  (PG)  cleave  glycosidic  bonds  with  the  aid 
of  water,  while  pectin  lyase  (PL)  and  pectate  lyase  (PAL) 
act  by  the  transelimination  mechanism.  Pectin  and  pectate 
lyases  are  also  called  pectin  and  pectate  transeliminases 
(PTE  and  PATE) (Whitaker,  1972). 

2.2.3  Mechanism  of  Action  of  Pectinases 

Generally  there  are  three  possible  ways  in  which 
enzymatic  degradation  of  alpha-1-4  linkages  can  occur  in 
pectin:  combined  action  of  pectinesterase  (PE)  and 
polygalacturonase  (PG) ;  hydrolytic  polymethylgalacturonase 
(PMG)  reaction;  and  pectin  transeliminase  (PTE)  reaction 
(Seegmiller  and  Jansen,  1952) .  In  the  process  of  enzyme 
action,  there  is  a  reduction  in  viscosity  depending  on  the 
temperature,  amount  of  enzyme  and  the  type  of  citrus  juice. 
Slowly  the  fine  haze  in  the  juice  begins  to  agglomerate  and 
form  a  floe  which  eventually  settles  out  leaving  the  clear 
supernatant  juice.  Often,  however,  even  the  clear 
supernatant  has  some  suspended  material  in  it.  Research  done 
by  Yamasaki  et  al.   (1967),  Ishii  and  Yokotsuka  (1971,  1972, 
1973)  and  Endo  (1964)  have  contributed  somewhat  to  the 
understanding  of  the  mechanism  of  action  of  pectinases. 
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Soluble  pectin  was  found  to  act  as  a  protective  colloid,  and 
partial  hydrolysis  of  the  pectin  would  cause  insoluble  and 
finely  divided  particles  to  flocculate.  Experiments  with  the 
haze  compounds  show  it  to  be  a  protein-carbohydrate  complex 
containing  about  36%  protein  (Yamasaki  et  al.,  1967).  At  a 
pH  of  3.5,  the  surfaces  of  these  protein  molecules  and 
pectin  as  well  as  other  carbohydrates  are  negatively 
charged.  Underneath  the  negative  charge  coat,  there  are 
postively  charged  proteins.  Partial  hydrolysis  of  this 
negative  coat  leads  to  the  exposure  of  positively  charged 
surfaces.  In  this  state,  the  system  has  unlike  charges 
attracting  one  another  and  a  floe  could  form. 

Cloud  from  orange  juice  contains  mainly  pectins, 
proteins  and  lipids  (Baker  and  Bruemmer,  1969) .  Hesperidin 
has  been  found  in  the  crystalized  form  in  the  orange  cloud 
(Mizrahi  and  Berk,  1970) .  Pectinesterase  is  the  main  enzyme 
which  determines  orange  cloud  flocculation  and  the  stability 
of  the  cloud  seems  to  be  a  function  of  the  ratio  between  the 
amount  of  flavonoid  hesperidin  and  the  molecular  weight  of 
the  polymer  pectin  (Ben-Shalom  et  al.,  1985). 

2.2.4  Composition  of  Citrus  Juices 

In  order  to  use  a  model  solution  which  would  simulate  a 
citrus  juice  (especially  orange  juice) ,  the  composition  of 
the  juice  has  to  be  known.  The  composition  of  citrus  juice 
has  been  studied  by  Cohen  et  al.   (1984),  Petrus  and 
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Vandercook  (1980).  Cohen  et  al.   (1984)  reported  an  acidity 
content  of  0.75  to  1.52  (%w/w)  for  Israeli  orange  juice  as 
citric  acid  anhydrous.  The  second  highest  acid  concentration 
was  that  of  ascorbic  acid  being  35  to  70  mg  vitamin  C/100  ml 
of  juice,  reported  by  Beacham  and  Bonney  (1937)   in  Florida. 
These  workers  reported  the  same  range  for  all  orange 
varieties  and  found  no  varietal  differences  in  vitamin  C 
levels.  However,  subsequent  studies  conducted  on  Florida 
orange  juices  by  Nagy  and  Smoot  (1977)  showed  that  juices 
from  early  season  (November  to  January,  mostly  Hamlin)  and 
mid  season  (January  to  March,  mostly  Pineapple)  contained 
more  vitamin  C  than  late  season  juices  (April  to  July, 
mostly  Valencia).  Birdsall  et  al.   (1961)  reported  that  juice 
from  California  naval  oranges  contained  about  15  mg/100  g 
more  total  vitamin  C  than  the  juice  from  Valencia  oranges. 

2.3  Modeling  of  pH  using  Reaction  Equilibria  in  The 
Simulation  Model  Systems 

2.3.1  Electrolyte  Solutions  and  Their  Modeling 

The  relative  importance  of  studies  of  electrolyte 
solutions  is  increasing  because  applications  to  phase 
equilibria  for  chemical  process  design  are  numerous  and 
industrially  important.  Process  design  requires 
applicability  in  a  large  range  of  concentrations, 
temperatures  and  pressures  for  complex  systems.  Maurer 
(1983)  gave  an  excellent  review  of  electrolyte  solutions. 
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In  the  case  of  aqueous  electrolytes,  the  solute  is 
polar  enough  to  allow  dissociation  of  electrolytes.  In  an 
ideal  situation,   it  may  be  assumed  that  a  particular 
electrolyte    completely  dissociates  in  the  liquid  phase. 
Deviations  from  ideality  are  mostly  derived  from  excess 
Gibbs  energy  at  given  temperature,  pressure  and  composition 
or  from  the  Helmholtz  energy  at  given  temperature,  volume 
and  composition  for  the  solution  of  "real"  species.  Choice 
of  the  reference  state  and  scale  of  concentration  are  often 
reviewed  and  depend  on  the  suitability  for  specific 
applications  (Renon,  1981;  Stokes,  1979).  The  properties 
which  can  be  derived  are  mainly  activity  coefficients, 
excess  partial  molar  volumes,  enthalpies,  entropies  and  heat 
capacities,  depending  on  the  availability  of  specific  data. 

For  calculation  of  the  mean  activity  coefficients, 
Meissner  and  Kusik  (1972),  Bromley  (1973)  and  Pitzer  (1973) 
have  proposed  several  models  based  on  observations  of 
experimental  behavior.  The  most  important  contribution  was 
that  of  Debye-Huckel  who  established  the  theoretical 
limiting  law.  Most  models  are  based  on  extended  Debye-Huckel 
equation  by  addition  of  a  polynomial  of  ionic  strength  or  by 
consideration  of  solvation  equilibria  in  order  to  extend 
their  applicability  to  higher  ionic  strength.  Single  ion 
activity  coefficients  can  be  calculated  by  the  Debye-Huckel 
expression  if  the  approximate  effective  radius  of  the 
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particular  ion  is  given  in  aqueous  solutions  (Kielland, 
1937)  . 

2.3.2  Equilibrium  Reactions  in  the  Liquid  Phase 

In  an  electrolyte  system  under  pressure  of  carbon 
dioxide,  a  number  of  major  reactions  including  dissociation 
reactions  occur  in  the  liquid  phase.  Considering  the 
composition  of  a  citrus  juice,  several  model  juice  systems 
can  be  used  to  simulate  a  citrus  juice.     Pure  water,  an 
aqueous  solution  of  citric  acid,  an  aqueous  solution  of 
ascorbic  acid  or  a  combination  of  all  these  can  be  used  with 
supercritical  carbon  dioxide.  The  dissolved  carbon  dioxide 
in  the  solution  and  the  presence  of  acids  change  the  pH  of 
the  model  systems.  The  significant  reactions  contributing  to 
the  lowering  of  the  pH  of  the  solution  are  the  dissociation 
of  ascorbic  and  citric  acids  and  the  formation  of  carbonic 
acid  as  a  result  of  carbon  dioxide  being  dissolved  in  the 
solution.  Figures  2-1  to  2-3  show  the  equilibrium  reactions 
occur ing  in  the  liquid  phase  of  the  citrus  model  systems 
considered. 

2.3.3  Effect  of  Dissolved  Carbon  Dioxide  on  pH  of  Pure  Water 
Although  the  dissociation  constant  data  of  carbonic 

acid,  produced  when  carbon  dioxide  is  dissolved  in  water,  is 
of  great  fundamental  importance,  from  a  practical  point  of 
view  the  true  acidity  expressed  in  terms  of  pH  is  more 
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useful.  Byke  (1932)  calculated  the  pH  values  for  solutions 
of  carbon  dioxide  in  pure  water.  For  a  partial  pressure  of 
zero  to  1.01x10  1  MPa  (atmospheric  pressure)  of  carbon 
dioxide.  In  his  calculations,  Byke  assumed  that  the 
dissociation  constant  remained  unchanged  with  process 
pressure.  However,  dissociation  constant  can  be  a  function 
of  temperature,  pressure,  and  composition.  Values  for  pH  of 
solutions  of  carbon  dioxide  in  water  at  pressures  above 
atmospheric  pressure  have  been  determined  colorimetrically 
by  Moore  and  Buchanan  (1930) .  These  measurements  were  made 
in  an  attempt  to  explain  the  action  of  carbon  dioxide  in 
carbonated  beverages  in  inhibiting  the  growth  of  bacteria. 
The  maximum  acidity  of  pure  water  was  reached  with  a  pH  of 
3.3  at  9.60xl0"1  MPa  when  the  temperature  was  25°C.  The 
conclusion  was  that  the  action  of  carbon  dioxide  solutions 
on  bacteria  was  due,  in  part  at  least,  to  some  factor  other 
than  the  lowering  of  the  solution  pH  (Moore  and  Buchanan, 
1930) . 

In  order  to  calculate  the  pH  more  accurately,  the 
variations  of  thermodynamic  properties  with  process  pressure 
have  to  be  taken  into  account,  namely,  the  variatons  of 
molar  volumes,  equilibrium  constants  and  activity 
coefficients. 
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2.3.4  Prediction  of  Thermodynamic  Properties  of  Aqueous 
Electrolytes  at  High  Pressures 

2.3.4.1  Pressure  dependance  of  molar  volume 

Upon  addition  of  an  electrolyte  to  water,  short-range 

order  among  the  solvent  molecules  is  destroyed  in  the 

vicinity  of  the  ions  by  local  collapse  of  the  water 

structure,  which  permits  H20  dipoles  to  orient  in  primary 

and  more  disordered  secondary  coordination  shells  about  the 

aqueous  species  (Bockris,  1949;  Samoilov,  1957;  Duncan  and 

Kepert,  1959) .  The  local  disruption  of  hydrogen  bonding 

among  the  solvent  dipoles  by  ion-solvent  interaction  is 

accompanied  by  a  decrease  in  entropy  as  well  as  volume. 

Relatively  few  equations  of  state  have  been  proposed  that 

would  account  for  the  volume  changes  accompanying  ion 

dissociation  and  ionization  as  a  function  of  pressure.  These 

equations  have  been  reviewed  by  Hamman  (1963)  and  more 

recently  by  Gancy  (1972),  Brummer  and  Gancy  (1972),  Helper 

and  Wooley  (1973),  Olofsson  and  Helper  (1975),  and  Hemmes 

(1972) . 

A  general  equation  of  state  for  aqueous  electrolytes 
was  described  by  Helgeson  and  Kirkham  (1976)  taking  into 
account  the  combined  effect  of  volume  loss  due  to  local 
collapse  of  the  solvent  structure  (Av°c)  and  that  relating 
to  the  solvation  process  (Av°s)  ;  that  is 


C  8 


(2-3) 
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Differentiation  of  the  above  equation  with  respect  to 
pressure  at  constant  temperature  leads  to 


{^tf~)t=  -k°=-w*  (2_4) 


where  k    represents  the  standard  partial  molar  isothermal 
compressibility,  and  <■>  and  N  are  the  Born  parameters  which 
can  be  determined  experimentally  knowing  the  partial  molar 
Gibbs  free  energy  of  solvation  for  a  particular  ion.  Using 
experimental  data,  the  isothermal  compressibility  of  an  ion 
can  be  calculated  from  the  above  equation  and  then  the 
change  in  the  molar  volume  can  be  obtained  by  integrating 
the  above  equation  from  the  standard  pressure  to  the 
operating  pressure.  Partial  molar  compressibilities  and 
partial  molar  volumes  are  given  for  a  number  of  uni-  and 
multivalent  ions  in  Helgeson  and  Kirkham  (1974a)  from  1  to  5 
kilobars  and  25  to  500  C. 

Knowing  the  partial  molar  volume  of  an  ionic  species, 
the  change  in  the  total  reaction  volume  associated  with  the 
dissociation  of  an  electrolyte  can  be  calculated  as  the 
difference  between  the  partial  molar  volumes  of  products 
minus  those  of  reactants: 
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change  of  molar  volume  of  reaction 
stoichiometric  coefficient  of  species  i 
partial  molar  volume  of  species  i 

The  changes  of  reaction  volume    are  given  for  a  number 
of  weak  acids  in  Owen  and  Brinkley  (1941) ,  Disteche  (1965) , 
Hamann  (1957) ,  Redlich  and  Bigeleisen  (1942) ,  Harned  and 
Owen  (1958)  and  Bodansky  and  Kauzmann  (1962) . 

2.3.4.2  Pressure  dependence  of  equilibrium  constant 

For  the  following  general  dissociation  reaction  of  an 
acid,  the  dissociation  equilibrium  is  characterized  by  the 
dissociation  constant,  K: 

HA    ■**   H*  +  A~ 

The  dissociation  constant  (also  refered  to  as  the 
equilibrium  constant)   is  the  quotient  of  the  activities,  af, 
of  the  chemical  species  involved  in  the  reaction.  The 
activity  is  a  measure  of  the  chemical  potential  of  a 
particular  substance  and  can  be  expressed  as  the  product  of 


where  AV°R  = 
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the  molar  concentration  and  the  corresponding  activity 
coefficient: 


K-U  a  - 


aH*  aA-  _  XH*  XA- 


(2-6) 


HA 


HA 


HA 


where  K„  =  equilibrium  constant  of  reaction  n 
a,-  =  activity  of  species  i 
Xj  =  molar  concentartion  of  species  i 
Y,  =  activity  coefficient  of  species  i 

The  equilibrium  constant  is  a  function  of  temperature, 
pressure  and  composition.  Although  the  equilibrium  constants 
for  a  number  of  dissociation  reactions  are  available  at 
mostly  the  standard  conditions  of  25°C  and  atmospheric 
pressure,  little  or  no  data  are  available  at  other 
temperatures  and  at  higher  pressures.  It  is  therefore 
necessary  to  develop  methods  of  determining  the  temperature 
and  pressure  dependence  of  the  equilibrium  constant  using 
the  available  thermodynamic  data  together  with  thermodynamic 
formulas. 

The  following  equation  is  the  mathematical  expression 
of  the  Le  Chatelier's  principle,  stating  the  pressure 
dependence  of  the  equilibrium  constant  at  constant 
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temperature (T)  and  constant  composition (x) ,   (Suzuki  and 
Taniguchi,  1972) : 


l^xT=^  (2-7) 
dP     x' T  RT 


where  K  =  equilibrium  constant 
R    =  ideal  gas  constant 
P    =  pressure 
T    =  temperature 
Av  =  change  of  reaction  volume 


2.3.4.3  Pressure  dependence  of  activity  coefficient 

Knowing  or  calculating  the  mean  activity  coefficients 
of  the  ionic  or  molecular  species  involved  in  the 
dissociation  reactions  at  standard  conditions,  the  activity 
coefficient  can  be  estimated  at  higher  pressures  at  constant 
temperature  and  composition  by  the  following  relation 
derived  from  Gibbs-Duhem  expression  (Disteche,  1972) : 
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dP 


X,  T 


RT 


(2-8) 


where  y=activity  coefficient 
P  =  pressure 

AV=  change  in  the  partial  moler  volume  of  the  ion 
R  =  ideal  gas  constant 
T  =  temperature 
x  =  composition. 

The  pH  of  an  electrolyte  solution  is  related  to  the 
activity  of  hydrogen  ions  in  the  solution.  The  equilibrium 
concentrations  of  ions  and  other  chemical  species  can  be 
calculated  through  expressions  that  relate  concentration  to 
the  thermodynamic  parameters  such  as  the  equilibrium 
constant,  the  compressibility  coefficient,  the  molar  volume 
and  the  activity  coefficient.  The  values  of  these  parameters 
at  higher  temperatures  and  pressures  can  be  estimated  by 
using  data  at  standard  conditions  given  in  the  literature 
along  with  the  appropriate  thermodynamic  relation.  The  next 
chapter  discusses  the  procedure  used  to  calculate  the 
equilibrium  concentrations  and  subsequently  the  pH  of  the 
electrolyte  solutions  used  as  the  citrus  simulation  systems 
in  this  study. 


CHAPTER  3 
MATERIALS  AND  METHODS 

This  chapter  gives  the  procedure  for  developing  the 
thermodynamic  models  for  the  four  simulation  systems  of  pure 
water-carbon  dioxide,  ascorbic  acid-water-carbon  dioxide, 
citric  acid-water-carbon  dioxide  and  ascorbic  acid-citric 
acid-water-carbon  dioxide.  In  developing  the  pH  model,  the 
model  parameters  such  as  the  activity  coefficients  and  the 
equilibrium  constants  were  first  calculated  and  later  used 
in  a  series  of  algebraic  equations  to  calculate  the 
equilibrium  concentrations  of  the  chemical  species  involved 
in  the  solutions.  The  pH  of  the  simulation  solutions  were 
later  calculated  from  the  equilibrium  concentration  of 
hydrogen  ions.     Once  the  pH  values  were  calculated  for  every 
simulation  system  at  various  process  conditions,  they  were 
compared  with  actual  experimental  pH  values  for  all  four 
systems  as  well  as  for  single  strength  orange  juice. 

3.1  Development  of  the  Thermodynamic  Models 
3.1.1  Calculation  of  Activity  Coefficients 
3.1.1.1  Pierotti  et  al.  method  for  carbonic  acid 

The  activity  coefficient  of  the  carbonic  acid  in  the 
mathematical  model  was  calculated  using  the  modified 
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equation  of  Pierotti  et  al.   (1959).     This  equation  relates 
the  infinite  dilution  activity  coefficients  to  the  molecular 
structures  of  the  solute  and  solvent  molecules  through  the 
number  of  carbon  atoms  for  both  species  as  shown  in  equation 
(3-1). 

logYVAi.a+B2-~+~+^(^-^)2+-^  (3-D 

2         1  2 


where  subscript  1  refers  to  solute  and  2  to  solvent  and 
Y*,        =        activity  coefficient  of  solute  (carbonic  acid) 
A1  2      =        coefficient  which  depends  on  nature  of  solute 

and  solvent  functional  groups 
B2         =        coefficient  which  depends  only  on  nature  of 

solvent  functional  groups 
=  coefficient  which  depends  only  on  solute 

functional  groups 
D         =         coefficient  independent  of  solute  and  solvent 

functional  groups 
F2       =         coefficient  which  essentially  depends  on  nature 

of  solvent  functional  groups 
N1,N2=  number  of  carbon  atoms  in  solute  and  solvent, 

respectively. 
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The  numerical  values  for  the  parameters  in  equation 
(3-1)  were  taken  from  Pierotti  et  al.   (1959)  parameter  table 
which  has  been  developed  for  a  number  of  organic  acids  and 
solvents.  Table  3-1  gives  the  parameter  values  for  the 
calculation  of  the  activity  coefficient  of  carbonic  acid  in 
water . 


Table  3-1.  Pierotti  parameters  for  calculating  the  infinite 
dilution  activity  of  carbonic  acid  in  water.  


A1,2 

-1.00 

0.662 

0.490 

D 

0.0 

F2 

0.0 

Hi 

1.0 

N2 

0.0 

3.1.1.2  Debye-Huckel  method  for  individual  ions 

Activity  coefficients  for  the  single  ions  in  the 
solution  were  calculated  using  an  extended  Debye-Huckel 
equation  (Pitzer,   1973) .  The  equation  is  based  on  the 
inter ionic  attraction  effects  of  the  aqueous  solution, 
expressed  as  the  ionic  strength,  /x,  of  the  solution. 


Equation  (3-2)  holds  for  the  range  of  ionic  strengths 
(0</x<0.7)  : 

0. 512  Z,ZBv^T 
-logy  .=  

l+aR/\i 


where  a  =  diameter  of  ion  in  Angstrom  (A) 

B  =  a  term  depending  on  temperature  and 

dielectric  constant  of  the  solvent 
Z-=  ionic  charge  of  species  i 
H  =  ionic  strength  of  solution 
The  approximate  values  for  the  effective  ionic  radii  of 
inorganic  and  organic  ions  are  given  in  the  following  table 
(Kielland,  1937) : 


Table  3-2.  Approximate  effective  ionic  radii  of  Inorganic 
and  organic  ions  


Ion  Species 

Approximate  Ionic  Radii 
(Angstroms) 

H+ 

9.0 

OH" 

3.5 

HCO/ 

4.0 

CO," 

4.5 

HAsc" 

3.5 

Asc 

4.5 

H,Cit" 

3.5 

HCit"" 

4.5 

Cit" 

5.0 
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3.1.1.3  UNIFAC  method  for  ascorbic  and  citric  acids 

The  UNIFAC  (UNIQUAC  Functional  Group  Activity 
Coefficients)  method  was  used  to  estimate  the  activity 
coefficients  of  ascorbic  and  cirtic  acids.  This  method  is  a 
group  contribution  method  that  uses  existing  phase 
equilibrium  data  and  applies  them  to  the  structural  groups 
of  various  compounds.  In  the  UNIFAC  method  the  logarithm  of 
the  activity  coefficient  is  estimated  to  be  the  sum  of  two 
contributions:  a  combinatorial  part,  essentially  due  to 
differences  in  size  and  shape  of  the  molecule  in  the 
mixture,  and  a  residual  part,  essentially  due  to  energy 
interactions  as  shown  in  equation  (3-3) : 

In  (yi)=ln(Yic)  +ln(Yi*)  (3-3) 


where      YjC=  combinatorial  activity  coefficient 
Y,R=  residual  activity  coefficient 
The    combinatorial    part    of    the    above    equation  was 
calculated  using  equations   (3-4)   to   (3-9)    from  Magnussen  et 
al.  (1981): 
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lny1c-ln-^+-|<7iln-3r*Ji-—  (Vi^)  (3"4) 

1  x1     2    1      (j^     1    xz     1  1     2  2 


and 


lny2c=ln^  +  4g2lnA+i    A(x_z  +x  j  )  (3~5) 
2  x2     2    2      <|>2     2    x2      1  1     2  2 


where 


4>1  = 

e1= 


<3Ti*i 


and 


 £2*2  


e2= 


qxx^q2x2 


(Noting  that  0,  +  <p2  =  1  and  9,  +  e2  =  1.) 


(3-6a,b) 


(3-7a,b) 
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V"§  (JV-Ok)- ;  z=10 


(3-8a,b) 


J2=f  (r2-g2)-(r2-l) 


(3-9a,b) 


where  xi'x2  =roole  fraction  of  component  1  and  2 

v  =number  of  groups  in  molecule  1  or  2 

Rk        =the  van  der  Waals  volume  for  group  k 

Qk        =the  van  der  Waals  surface  area  for  group  k 

k         =group  number  (assigned)  from  table  of 

group  volumes  and  surface  areas  (Table  3-3)  . 

The  group  volumes  and  surface  areas  for  the  combinatorial 
activity  coefficients  were  taken  from  the  UNIFAC  tables  of 
vapor-liquid  equilibrium  data  by  Magnussen  et  al. (1981)  and 
are  given  in  Table  3-3. 
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Table  3-3. Group  volume  and  surface-area  parameters  for  citric 
and  ascorbic  acids  used  in  the  UNIFAC  method  


Acid/Parameter 

k 

Qk 

Citric (-COOH) 

43 

1.3013 

1.224 

Citric  (-CH2) 

2 

0.6744 

0.540 

Citric ( -OH) 

15 

1.000 

1.200 

Citric (-C-\ 

4 

0. 2195 

0 

A  c?  r*OT*Vi  "i  c  (  — OT-T  ^ 

1 5 

1 . 000 

1.200 

Ascorbic 

(-C=C-) 

9 

0.6605 

0.485 

Ascorbic ( -C-H) 

3 

0.4469 

0.228 

Ascorbic 
J  (-CCOO-) 

23 

1.6764 

1.420 

The  residual  part  of  the  equation  for  calculating  the 
activity  coefficients  of  ascorbic  and  citric  acids  was 
calculated  using  equations  (3-10)  to  (3-12) : 


lnYl*=y  v  <1)  [lnlVlnTic'l)  ] 

^  (3-10a,b) 

lnY2*=£v k{2)  [lnr^-lnTic'2)  ] 
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where 


inrJt-oJt[i-in(j;e.*Jlk)-2  (0,#jj;e^j]  o-n) 


and  m  and  n=l, 2 , . . . ,N(all  groups) 


•*  Ew.    '     *  '  ' 

j=l,2,...,M;  n=l,2,...,N,  Also, 

j        =    component  number  (1  or  2  for  a  binary  mixture) 
M       =    total  number  of  components (=2  for  a  binary 

mixture) 
n        =    group  number 
N        =    total  groups 
6m      =    group  surface  area  fraction 
Xra      =    group  fraction (For  example,  for  a  equimolar 

mixture  of  butane  and  hexane,  the  group 

fractions  for  CH3  and  CH2  are  4/10  and 

6/10,  respectively.) 
Xj        =  mole  fraction  of  molecule  j  in  the  mixture 
v(J)m  =  number  of  groups  of  type  m  in  molecule  j 
^      =  expt-a^/T) 
T  =  temperature (K) 
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The  parameter  a^  is  a  group  interaction  parameter  which 
is  a  measure  of  the  difference  in  interaction  energy  between 
groups  n  and  m  given  in  Magnussen  et  al.   (1981) . 

The  solubility  of  carbon  dioxide  in  the  solution  at  the 
process  temperature  and  pressure  was  approximated  using  vapor- 
liquid  equilibrium  data  for  carbon  dioxide  and  water. 
Solubility  values  for  carbon  dioxide  has  been  published  by 
Wiebe  and  Gaddy  (1940)  for  the  C02-H20  binary  system  between 
12  and  100°C  and  at  pressures  from  2.53  to  70.70  MPa.  The 
solubility  of  carbon  dioxide  at  atmospheric  pressure  was 
calculated  using  equation  (3-13) ,   (Wiebe  and  Gaddy,  1940) : 

J?  InS  =  A  +  B/T  +  C  InT  +  D  T  (3-13) 


where  S=mole  fraction  solubility  of  carbon  dioxide  in  pure 
water  at  one  atmosphere  (1.0136xl03  KPa) 
A=-1330.03  J/ (K. mole) 
B=7.2733  KJ/mole 
C=180.295  J/ (K. mole) 
D=-0. 009174  J/ (K2. mole) 
R=8. 31957  J/ (K. mole) 

T=absolute  temperature  (273  K<  T  <  353  K) 


An  interpolation  program  was  written  to  calculate  the 
solubility  at  the  exact  temperature  and  pressure  of  the 
supercritical  process.  The  parameters  of  the  thermodynamic 
model  were  either  taken  from  data  available  in  thermodynamic 
literature  or  calculated  using  the  reaction  equilibria 
formulae.  These  parameters  included  the  equilibrium  constants, 
k,  the  molecular  and  single  ion  activity  coefficients,  y, 
compressibility  coefficients,  k,  and  changes  of  reaction 
volumes,  AV. 

3.1.2  Calculation  of  pH  of  the  Simulation  Systems 

The  pH  of  the  model  solutions  was  calculated  using 
equation  (3-14)  expressing  the  pH  as  the  negative  logarithm  of 
the  activity  of  the  hydrogen  ions: 


pH=  -loga^  =  -logxH+-logy 


B* 


(3-14) 


where  aH+=activity  of  hydrogen  ions 

xH+=concentration  of  hydrogen  ions 
yH+=activity  coefficient  of  hydrogen  ions. 
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3.1.3  Systematic  Method  for  Calculating  Equilibrium 
Concentrations 

A  systematic  method  was  used  to  determine  the  equilibrium 
concentrations  of  the  unknown  species  in  the  model  solution. 
The  first  step  is  to  define  a  set  of  equations  for  all  the 
reactions  that  appear  to  have  a  major  effect  on  the 
concentration  of  hydrogen  ions.  These  equations  included 
expressions  for  the  equilibrium  constant  in  the  major 
equilibria  for  the  particular  model  system,  the  mass-balance 
equation,  and  the  charge-balance  equation.  Mass  balance 
equations  were  algebraic  expressions  relating  the  equilibrium 
concentrations  of  the  various  species  to  one  another  and  to 
the  analytical  concentrations  of  the  substances  present  in  the 
solution.  The  charge-balance  equation  expresses  the 
relationship  that  the  concentrations  of  the  cations  and  anions 
must  be  such  that  the  solution  is  electrically  neutral.  All 
of  these  expressions  will  result  in  a  number  of  independent 
equations  which  will  be  equal  to  the  number  of  unknown 
concentrations . 

3.1.3.1  Ionization  constant  of  water 

Marshall  and  Frank  (1981)  gave  the  ion  product  of  water 
for  a  range  of  temperatures,  shown  in  Table  3-4: 
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Table  3-4.  Ionization  constant  for  water  (KJ 


-log 

Temperature ,  °C 

-log  K, 

Temperature ,  °C 

14.9435 

0 

13.8330 

30 

14.7338 

5 

13.6801 

35 

14.5346 

10 

13.5348 

40 

14.3463 

15 

13.3960 

45 

14.1669 

20 

13.2617 

50 

14.0000 

24 

13.1396 

55 

13.9965 

25 

13.0171 

60 

3.1.3.2  Numerical  values  of  equilibrium  constants 

The  equilibrium  (dissociation)  constants  of  all  of  the 
reactions  involved  in  the  model  solutions  were  found  in  the 
literature  and  are  given  in  Table  3-5  (Meites,  1963) : 
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Table  3-5.  Dissociation  constants  for  carbonic,  ascorbic  and 
citric  acids  


Acid 

Dissociation  Constant 

Carbonic  (HC03" ) 

K^s.eixio"11 

Carbonic (H2C03  to  C03") 

K2=4.31xl0"7 

Carbonic  (H2C03  to  C02) 

K3=4.42xl0"7 

Ascorbic  (C6H806) 

K4=7.94xl0"5 

Ascorbic  (C6H706  ) 

K5=l.  62X10"12 

Citric  (C6H807) 

1^-7.45x10"* 

Citric  (C6H707  ) 

K^l.TSxlO5 

Citric  (C6H607") 

Kg=4.02xl0"7  J 

3.1.3.3  Set  of  equations  for  the  pure  water-carbon  dioxide 
system 

The  chemical  reactions  involved  in  the  binary  solution 
of  pure  water  and  carbon  dioxide  involved  the  dissociation  of 
water  and  the  three  dissociation  reactions  of  the  carbonic 
acid.  These  four  equations  were  derived  using  the  expressions 
for  the  equilibrium  constants  shown  in  equations  (3-15a)  to 
(3-15d) : 


K  _  aWaOH- 
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(3-15a) 


_  acq,2~  •  aw* 


(3-15b) 


awoof  •  a/r 

aH2COi 


(3-15C) 


if3  =  a^°,a°^  (3-15d) 

aH2CX53 


In  the  above  equations,  the  activities  can  be  written  as 
products  of  molar  concentrations  and  activity  coefficients,  as 
shown  in  equations  (3-16a)  to  (3-16d) : 


W H*  •  Y OH')  XB*  •  XOH'  ~  Kv  ~  0 


(3-16a) 
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(3-16b) 


(Y^-Yac^Ox^.x^-  -  (Y^co.-JC,)^^  =  0  (3-16c) 

(Ycd,)^  -  (Y^cDi-^J^co,  ■  0  (3-16d) 

The  two  other  equations  that  relate  the  concentrations  of 
different  species  in  the  solution  together  are  the  mass 
balance  and  the  charge  balance  equations  (3-17)  and  (3-18) : 


Carbon  mass  balance  equation: 


xco2  +  *h2co3  +  xhco;  +  xcol'  ~  sco2  =  0  (3~17) 


Charge  balance  (electroneutrality)  equation: 


xr  ~  xoh-  ~  xhco;  ~  2xcol~  ■  0  (3-18) 
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3.1.3.4  Set  of  equations  for  the  ascorbic  acid-water-carbon 
dioxide  system 

In  addition  to  the  equilibrium  constant  equations  in  the 

pure  water-carbon  dioxide  system,  two  more  equations  can  be 

written   for  the  Ascorbic   acid-water-carbon   dioxide  system 

which  take  into  account  the  dissociation  reactions  of  ascorbic 

acid.  The  equations  for  charge  and  mass  balances  have  to  be 

modified  because  of  additional  ions  present  in  the  solution. 

The  equilibrium  constant  equations  are  given  below: 


*4  = 


aH*  •  ac6n,o; 


(3-19) 


(3-20) 


Substituting  the  activities  with  the  product  of  molar 
concentrations  and  the  activity  coefficients,  equations  (3-21) 
and  (3-22)  were  obtained: 

WJT*'  Y csh,o^  xn*  •  xcsh,o;  ~  KiY csh„o6  •  XC6HaOs  =  ®  (3~21) 
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( Yjr  •  Ycs^o« " )  •  xh*  •  xcs/^o|~      %s  •  Ycffaol  •  xc6HjC£  ~  ®  (3-22) 


The  mass  balance  equation  (3-23) : 

XC02  +  XH2COi  +  -^ifOOj"  +  XCC>1~  ~  ^XC6Hg06 


The  charge  balance  equation  (3-24) : 


XM*  ~  xoh-  ~  xhcoI  ~  xctn,ol  ~  2xc{Hsot  ~  2xcol'  =  0  (3  "2  4) 


3.1.3.5  Set  of  equations  for  the  citric  acid-water-carbon 
dioxide  system 

In  this  ternary  system,  there  will  be  three  other 
equilibrium  constant  equations,  as  well  as,  different  mass  and 
charge-balance  equations  due  to  new  chemical  species  in  the 
model  solution: 
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Equilibrium  constant  equations  (3-25)  to  (3-27) : 


K6  =    a"* '  fa&g  (3-25) 


K,  -   ^V^'  (3-26) 


Ke  =       Z    ™  (3-27) 


By  substitutinq  the  activities  with  molar  concentrations 
and  the  activity  coefficients,  expressions  (3-28)  to  (3-30) 
were  obtained: 

CsH,0;  •  Y  H*)  XH*  •  XCtfyO;  "  K6  •  YCgHgO,  •  xc6h9o,  =  0  (3-28) 
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The  mass-balance  equation  (3-31) : 

XC02  +  XH2C03  +  XHCOi  +  ■XCO|"  "  ^XCsK807   +   ^^Cgi^O,"  + 
6XC6H60?-   +  6XC6HsO?—  5C0j   =  0 

The  charge-balance  equation  (3-32) : 

XH*~X0H-~XHC0;~XCsH,0;~2xCsHs0?-~3XCiH50}~~2XC0Z~  =   0  (3~32) 


3.1.3.6  Set  of  equations  for  the  ascorbic  acid-citric 
acid-vater-carbon  dioxide  system 

The    set    of    equations    resulting    from    the  equilibrium 

expressions  for  the  three  previous  simulation  systems  were  all 

used  in  developing  the  pH  model  for  the  quaternary  system  of 

ascorbic  acid-citric  acid-water-carbon  dioxide.  However,  the 

mass  balance  and  charge  balance  equations  would  be  different 


(3-29) 


(3-30) 
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since  both  ascorbic  acid  and  citric  acid  are  present  in  the 
solution: 


The  mass  balance  equation  (3-33) : 


x 


CO, 


+  x 


H-.CO 


j   +  XHCOi 


+  xcot 


~  6XC6Hi0f  + 


6XCSHSC^- 


~  6xcsHta,  +  6xctn,oJ-  + 


(3-33) 


+  6xc6«5o7- 


-  ^  =  0 


The  charge-balance  equation  (3-34) : 


XH*    XOH-    XHCCq    XCtH1Ol  ^■XCsHfol~ 


-x, 


-2x. 


2--3x. 


cshso; 


—2x. 


CO, 


=  0 


(3-34) 


In  developing  the  thermodynamic  models,  the  set  of  equations 
resulting  from  the  expressions  for  dissociation  constants  and 
mass  and  charge  balances  were  solved  simultaneously  for  the 
activity  of  hydrogen  ions  (H+)  using  the  Newton-Raphson 
iteration  method.  The  molar  concentrations  of  all  the  ionic 
and  non-ionic  species  in  the  solution  made  up  a  matrix  of 
eleven  elements.  By  assigning  initial  concentration  values  and 
inverting  the  coefficient  matrix,  the  iterated  concentrations 
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converged  to  a  set  of  new  values  corresponding  to  the 
equilibrium  concentrations  in  the  solution.  A  FORTRAN  program 
was  written  to  do  the  iteration  procedure  and  calculate  the 
concentration  of  the  hydrogen  ions  at  the  process  temperature 
and  pressure.  A  listing  of  this  program  is  given  in  the 
Appendix. 

3 . 2  Experimental 
A  series  of  experiments  were  performed  under  different 
process  conditions  with  different  levels  of  acid  concentration 
to  compare  the  pH  values  predicted  by  the  mathematical  models 
to  the  actual  pH  of  the  solution.  The  pH  values  were  measured 
at  high  pressures,  and  for  this  purpose  a  specially  designed 
pH  probe  was  used  which  could  withstand  high  pressures. 

3.2.1  High  Pressure  pH  Probe 

The  activity  of  hydrogen  ions  in  the  simulation  solutions 
was  measured  with  a  pH  measuring  system  consisting  of  a  glass 
electrode,  a  reference  electrode,  and  a  pH  meter.  The  reading 
on  the  pH  meter  is  the  glass  electrode  potential  which  is 
proportional  to  the  hydrogen  ion  activity  corrected  by  the 
reference  electrode  potential.  The  major  components  of  the 
glass  electrode  are  the  pH-sensitive  glass  bulb,  the  probe 
shell,  and  the  filling  solution.  The  glass  bulb  is  a  membrane 
which  shows  an  electric  potential  when  it  comes  in  contact 
with  hydrogen   ions.    Ordinary  pH  probes   have  very  limited 
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pressure  ranges  and  can  not  be  used  for  measurements  at  high 
pressures.  At  high  pressures  the  glass  bulb  and  the  probe 
casing  have  to  be  made  of  resistant  polymeric  material  to 
avoid  destruction  of  the  probe,  since  the  internal  pressure  of 
the  probe  would  be  lower  than  the  external  pressure  and  this 
leads  to  the  collapse  of  the  probe  casing.  The  high  pressure 
pH  probe  used  in  the  experiments  was  a  Markson  (Phoenix,  AZ.) 
with  a  gel  filled  combination  Ag/AgCl  electrode,  a  Teflon 
reference  junction,  and  a  stainless  steel  bulkhead  fitting. 
This  special  probe  was  originally  designed  for  taking  pH 
measurements  at  deep  places  of  oceans  where  the  probe  will  be 
under  high  pressures.  The  pH  meter  was  an  Orion  EA  920 
(Cambridge,  MA.). 

3.2.2  Experimental  Design 

The  experimental  design  chosen  to  determine  the 
combinations  of  model  variables  was  the  Box  and  Behnken(1960) 
design.  In  this  design,  a  fraction  of  the  three-level 
factorial  arrangements  was  used.  The  levels  of  each  of  the 
quantitative  factors  were  assumed  equally  spaced  and  they  were 
coded  -1,0,  and  1.  The  most  commonly  used  coding  scheme  was 
used  to  define  the  coded  variables,  xf,  which  were  defined  in 
standardized  form  as  shown  in  equation  (3-35) : 


55 

where  X,-  m  is  the  mean  of  the  Xuj  values  (u=l,  2 , . . .  ,N)  and  Si  is 
the  scale  factor. 

The  coded  variables  were  applied  to  the  experimental 
variables  which  included  six  levels  of  process  pressure  and 
three  levels  of  concentration  for  ascorbic  and  citric  acids 
each  and  for  every  simulation  system. 

3.2.3  Preliminary  Experiments 

A  set  of  preliminary  experiments  were  performed  to 
determine  the  effect  of  sugar  (sucrose)  on  the  solubility  of 
carbon  dioxide  and  the  subsequent  change  in  pH  of  the 
solution.  Other  experiments  were  performed  to  determine  the 
effect  of  process  temperature  on  the  pH  of  the  model 
solutions.  The  sucrose  concentrations  were  at  0  and  20%,  and 
the  process  temperatures  were  32,  37,  and  42°C.  These 
experiments  were  performed  for  six  levels  of  process  pressure, 
namely,   0,   0.690,   1.379,   2.758,   4.137,   5.861  MPa. 

3.2.4  Experiments  on  the  Effect  of  Pressure  and  Acid 
Concentrations  on  the  pH  of  Model  Solutions 

To  confirm  the  pH  values  predicted  by  the  thermodynamic 

model,  a  set  of  experiments  were  performed  with  a  solution  of 

citric  acid  (  certified  A.C.S.  reagent,  Registry  No.  876435 

provided  by  Fischer  Scientific) and  ascorbic  acid(  certified 

A.C.S.    Reagent,    Registry   No.    03912LX   provided   by  Aldrich 

Chemical    Company,    Inc.,    Milwaukee,    Wisconsin)    to  simulate 

citrus  juice.  The  concentration  levels  of  the  two  acids  and 
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process  temperature  and  pressure  were  varied  according  to  the 
Box-Behnken  design (  1960).  Table  3-6  shows  the  design 
variables  and  their  levels  in  the  experiments .  Concentration  of 
citric  acid  varied  from  0  to  1.17  grams  per  liter  of  solution 
and  that  for  ascorbic  acid  was  from  0  to  19.5  grams  per  liter 
of  solution.  Process  pressure  ranged  from    0  to  34.48  MPa. 


Table  3-6.  Coded  design  variable  settings  for  the  pH 
experiments  


Pressure, atm. 

Citric  Acid 
Concentration , g/L 

Ascorbic  Acid 
Concentration , g/L 

+1 

+1 

0 

+1 

-1 

0 

-1 

+1 

0 

 -  

-1 

0 

+1 

0 

+1 

+1 

0 

-1 

0 

+1 

-1 

0 

-1 

0 

+1 

+1 

0 

-1 

+1 

0 

-1 

-1 

0 

0 

0 
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In  the  experiments  with  pure  water,  a  solution  of  pH 
adjuster (  pHiX  Adjuster  provided  by  Orion  Research  Inc., 
Cambridge,  MA.)  was  used  to  stabilize  the  pH  readings,  since 
the  concentration  of  electrolytes  in  the  solution  would  be 
low.  A  high  pressure  vessel  was  used  as  a  closed  system  in  the 
high  pressure  experiments,  and  the  pH  was  monitored  by  the 
Orion  pH  meter  during  the  process.  Experiments  were  conducted 
in  a  constant-temperature  water  bath  (Blue  M,  Blue  Island, 
IL.).  Equilibration  time  was  a  minimum  of  2  hours,  determined 
in  the  preliminary  experiments.  Figure  3-1  shows  the  schematic 
of  the  experimental  set-up  for  measuring  the  pH  in  the 
supercritical  vessel.  To  take  the  pH  measurements,  40  ml  of 
solution  was  put  into  the  stainless  steel  pressure  vessel  and 
then  the  pH  electrode,  previously  fixed  to  the  lid  of  the  of 
vessel,  would  fit  into  the  head  of  the  vessel.  The  rubber  0- 
rings  on  the  lid  prevented  any  leaks  from  the  pressure  vessel 
and  before  each  run  the  whole  system  including  pressure 
valves,  pressure  hoses  and  the  pressure  vessel  were  purged  of 
air  by  C02  gas  to  minimize  the  oxidation  of  the  acid 
solutions.  The  pressure  of  the  vessel  was  regulated  by  a 
sensitive  pressure  valve  and  read  on  a  pressure  gage. 
Temperature  inside  the  vessel  was  maintained  by  the  water  bath 
and  the  readings  were  taken  by  a  mercury  thermometer  as  well 
as  the  pH  meter  thermocouple.  The  pH  electrode  was  kept  in  an 
Orion  storage  solution  to  prevent  drying  and  calibrated  by 
standard  buffer  solutions  before  each  experiment.    Once  the 
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solution  was  placed  in  the  pressure  vessel  and  subsequently 
into  the  water  bath,  15  minutes  was  allowed  to  pass  for  the 
solution  to  come  to  the  same  temperature  as  the  water  bath 
before  pressurization. 

3.3  Statistical  and  Sensitivity  Analyses 
The  pH  values  obtained  from  the  experiments  with  the 
effect  of  sugar  concentration  and  process  temperature,  were 
compared  using  a  t-test  statistics  with  a  0.5%  level  of 
significance.  The  null  hypothesis  was  that  the  difference  of 
mean  pH  values  for  each  set  of  treatments  was  not  different 
(d=  /i,  -  /i2  =  0)  .  The  test  statistic  was  determined  using 
equation  (3-36) : 


t  -  (3-36, 


where  davg  =  average  difference  of  pH  values  of  the  two 
treatments 

Ma      =  Mi  -  M2  =  difference  between  mean  pH  values 
sd      =  standard  deviation  of  paired  differences 
n        =  number  of  paired  samples 

The   correlation   coefficients    for   the   pH   models  were 
calculated  using  the  following  relation  (3-37) : 
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E^-P*)2  (3-37) 
£  (pH  -  pH)  2 

where  r  =  coefficient  of  correlation 

S(pH  -  pH)    =  regression  sum  of  squares 
S(pH  -  pH)2  =  total  sum  of  squares 
pH=  mean  pH  value 

A 

pH=  estimated  pH  value 
pH=  measured  pH  value 

Since  the  parameters  of  the  pH  models  developed  for  the 
four  simulation  systems  were  either  calculated  or  taken  from 
the  literature,  and  some  assumptions  were  made  in  estimating 
them,  their  effect  on  the  predicted  pH  values  were  studied. 
The  calculated  activity  coefficients  and  equilibrium  constants 
used  in  the  pH  model  were  decreased  or  increased  by  a  factor 
of  10%  to  determine  how  sensitive  the  model  is  to  variation  of 
these  parameters.  Also,  the  effects  of  carbon  dioxide 
solubility  values  as  well  as  the  influence  of  the 
concentrations  of  ascorbic  and  citric  acid  on  the  accuracy  of 
the  pH  models  were  also  studied. 


CHAPTER  4 
RESULTS  AND  DISCUSSION 

For  the  four  simulation  systems  of  pure  water-carbon 
dioxide,  ascorbic  acid-water-carbon  dioxide,  citric  acid- 
water-carbon  dioxide,  and  ascorbic  acid-citric  acid-water- 
carbon  dioxide,  the  activity  coefficients  and  the  equilibrium 
constants  were  first  calculated  and  later  used  in  the 
iteration  program  to  calculate  the  pH.  Whenever  there  were 
available  thermodynamic  data,  the  variations  of  the  activity 
coefficients  and  the  equilibrium  constants  with  pressure  were 
also  calculated.  The  values  of  pH  calculated  by  the 
thermodynamic  model  were  later  compared  with  the  pH  values 
measured  experimentally  for  each  simulation  model. 

4.1  Pure  Water-Carbon  Dioxide  Simulation  System 
4.1.1  Calculated  Activity  Coefficients 

The  activity  coefficients  for  the  ions  and  other  chemical 
species  in  the  water-carbon  dioxide  binary  system  were 
calculated  using  the  methods  discussed  in  the  previous  chapter 
and  are  given  in  Table  4-1  with  the  corresponding  method  of 
calculation. 
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4.1.2  Variation  of  Activity  Coefficients  with  Pressure 

Equation  2-4  in  chapter  2  relates  compressibility 
coefficient  to  change  in  partial  molar  volume  as  pressure  is 
increased  at  a  constant  temperature.  By  integrating  equation 
2-4  with  respect  to  pressure  and  using  the  data  given  in  the 
literature  for  compressibility  coefficients  and  partial  molar 
volumes  at  standard  conditions,  the  change  in  molar  volume  was 
calculated.  Assuming  compressibility  coefficients  stay 
constant  with  temperature  and  pressure,  the  partial  molar 
volume  of  the  ions  would  be: 


V±  =  10 


logy,*  ♦  *rP0-*rP 

*  1  2.303  (4-1) 


where    V^partial  molar  volume  of  species  i  at  pressure  P 
and  temperature  T 
Vj°=partial  molar  volume  of  species  i  at  standard 
pressure 

(101.36  KPa)  and  temperature  (25°  C) 
kT=isothermal  compressibility  coefficient 
P0=standard  pressure 
P  =process  pressure 
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The  change  in  the  molar  volume  of  the  ions  with  respect 
to  process  pressure  (AP  =  P  -  P0)  would  then  be  AV{  -  Vf  -  V,-0 

Table  4-2  gives  the  values  of  the  isothermal 
compressibility  coefficients  for  the  four  ions  present  in  the 
water-carbon  dioxide  system. 

Since  these  values  for  the  compressibility  coefficients 
were  given  in  cm3/mole.KPa,  integration  of  equation  2-4 
without  the  pressure  term  in  the  denominator  would  result  in 
the  following  equation: 


Using  the  data  given  in  Table  4-2,  the  change  of  molar 
volume  with  process  pressure  for  the  four  ions  in  the  solution 
are  given  in  the  following  equations: 


Vi  =  Vf  +  -  ic^P 


(4-2) 


aVh.  =  0  (cm*/mole) 


(4-3a) 


AVnu-  =  -2.9x10'*?  +  5.0xl0-2 


(cm2 /mole) 


(4-3b) 


nVHCO-  =  -4.Oxl0-5P  +  S.OxlO-3     (cm3  /mole) 


(4-3c) 
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Table  4-1.  Calculated  activity  coefficients  for  water-carbon 
dioxide  system 


tneniicai  bpecies 

taicuidLeu 
Activity 
Coefficient,  y 

MofhnH  nf 
rlB  L.I1CJU.  OX 

Calculation 

H+ 

0.830 

Debye-Huckel 

OH" 

0.760 

Debye-Huckel 

co3" 

0.360 

Debye-Huckel 

HC03" 

0.775 

Debye-Huckel 

H2C03 

1 

0.412 

Pierotti  et  al. 

Table  4-2.  Compressibility  coefficients  of  the  ions  in  the 
water-carbon  dioxide  system  


Ionic  Species 

Isothermal  Compressibility 
Coefficient, 
kT  ( cm3/mole  ( KPa ) 

H+ 

0 

OH" 

-2.94X10"1 

HC03" 

-4.05X10"2 

co3 " 

-7.60X10"1 

(Source:  Helgeson  and  Kirkham,  1976) 


AVm-  =  -7  .5x10-* P  +  7.5xl0-3     (cm3 /mole) 
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(4-3d) 


The  change  of  molar  volume  for  carbonic  acid  (AVH2C03)  is 
given  by  Disteche  (1975)  to  be  -27.0  cm3/mole  at  standard 
conditions.  Using  the  above  equations  and  the  change  of  molar 
volume  for  carbonic  acid,  the  variations  in  the  activity 
coefficients  with  pressure  in  the  water-carbon  dioxide  system 
were  formulated  using  equation  2-8  and  are  given  below: 


logy*.  =  -9.69x20"2  (4-4a) 

logYcy,-  =  -1.193X10-1  +  l.O3Oxl0_5P  -  2.575X10-*P2  (4-4b) 

logy =  -9  .69x1 0"2  (4-4c) 

logYtf2co3  "  -3.774X10-1  -  4.79x20-5P  (4-4d) 

logy^-  =  -1.1O7X20"1  +  1.421x20-6P  -  3.55Oxl0'9P2  (4-4e) 

logYc^--  =  -4.436X10"1  +  1.332x20'6P  -  6.661x20"8P2  (4-4f) 
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4.1.3  Variation  of  Equilibrium  Constants  with  Pressure 

The  change  in  the  reaction  volume  for  the  four  main 
equilibrium  reactions  in  the  water-carbon  dioxide  system  were 
found  in  the  literature  and  are  given  below:  (Meites,  1963; 
Disteche,  1972) 


^  =  - 


20 .3cm3 /mole 
27 .  0cm3 /mole 
29 . 0cm3  /mole 
27  .  Bern3 /mole 


Using  the  above  data  for  the  change  in  reaction  volume 
together  with  equation  2-7,  the  following  expressions  were 
obtained  for  the  equilibrium  constants  as  a  function  of 
process  pressure.  Equation  2-7  was  integrated  with  respect  to 
pressure  at  constant  temperature: 

logKw  =  -13.677  +  3.6xlO~5P  (4-5a) 
logiq  =  -10.252  +  5.0xlO-5P  (4-5b) 


logiC,  =  -6.366  +  5.0xlO~5P 


(4-5c) 
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logic,  =  -6.355  +  4.8xlO"5P 

The  above  expressions  for  activity  coefficients  and 
equilibrium  constants  as  functions  of  pressure  were 
substituted  in  the  set  of  simultaneous  equations  for  the 
water-carbon  dioxide  system  given  in  chapter  3.  Solution  of 
these  equations  by  the  Newton-Raphson  iteration  method  gave 
the  equilibrium  concentration  for  the  hydrogen  ions  in  the 
solution  from  which  the  pH  of  the  system  was  calculated. 

4.1.4  Comparison  of  the  Experimental  and  Calculated  pH  values 
4.1.4.1  Effect  of  extreme  process  pressures 

The  pH  values  up  to  a  pressure  of  34.48  MPa  predicted  by 
the  mathematical  model  are  compared  with  the  experimental  pH 
values  measured  by  the  high  pressure  probe  in  Table  A-l  (in 
Appendix  A)  .  Figure  4-1  is  the  plot  of  the  predicted  and 
measured  pH  values  at  different  process  pressures.  As 
indicated  by  this  plot,  the  measured  pH  values  are  closely 
predicted  by  the  pH  values  calculated  by  the  mathematical 
model  for  the  water-carbon  dioxide  system. 

The  greatest  lowering  of  pH  of  the  water-carbon  dioxide 
system  happened  at  pressures  up  to  about  5.516  MPa.  At  higher 
pressures  the  pH  of  the  solution  did  not  change  significantly 
with  further  pressure  increase.  This  trend  of  the  pH  of  the 
binary  system  was  observed  both  in  the  experimental  pH  values 
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as  well  as  in  the  pH  values  predicted  by  the  thermodynamic 
model.  From  a  pressure  of  0  to  5.516  MPa  the  variation  in  the 
measured  pH  was  a  drop  of  2.53,  whereas,  from  5.516  MPa  to 
34.48  MPa  the  pH  was  lowered  by  only  0.14. 

Using  color imetric  methods,  Moore  and  Buchanan  (1930) 
found  the  same  trend  of  pH  lowering  for  aqueous  solutions  of 
carbon  dioxide  upto  process  pressures  of  3.363  MPa,  and  at 
this  process  pressure,  the  measured  pH  of  the  aqueous  solution 
was  3.3.  Figure  4-2  compares  the  predicted  pH  values  for  the 
pure  water-carbon  dioxide  system  with  the  measured  pH  values 
in  this  study  as  well  as  the  measured  pH  values  measured 
colorimetrically  by  Moore  and  Buchanan  (1930) .  There  is  a  good 
agreement  between  the  predicted  and  the  measured  pH  values  by 
both  methods. 

The  sharp  decrease  followed  by  leveling  in  the  pH  values 
in  the  pure  water-carbon  dioxide  system  may  be  explained  by 
the  fact  that  since  the  pressure  vessel  is  a  closed  system, 
the  dissolution  of  carbon  dioxide  in  water  may  be  approaching 
a  saturation  limit  at  these  pressures. 

4.1.4.2  Effect  of  sugar  on  solubility  of  carbon  dioxide 

The  measured  pH  values  for  pure  water  and  an  aqueous 
solution  of  20%  sucrose  solution  are  qiven  in  Table  4-3. 
These  values  do  not  show  any  significant  difference  in  the  pH 
of  the  solution  implying  that  up  to  the  pressure  of  5.86  MPa, 
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sucrose  concentrations  of  up  to  20%  do  not  significantly 
affect  the  solubility  of  carbon  dioxide  in  the  solution. 

4.1.4.3  Effect  of  process  temperature  on  pH 

Table  4-4  gives  the  measured  pH  values  for  the  pure 
water-carbon  dioxide  system  at  three  process  temperatures  of 
32,  37  and  42°C  over  a  pressure  range  of  0  to  5.861  MPa.  Over 
this  experimental  temperature  and  pressure  range,  no 
significant  change  in  the  pH  values  was  observed.  However, 
the  slightly  lower  pH  values  at  the  lower  temperatures  may  be 
explained  by  the  fact  that  the  solubility  of  carbon  dioxide  in 
water  is  greater  at  low  temperatures  than  at  high  temperatures 
(Weibe  and  Gaddy,  1940) . 

4.2  Ascorbic  Acid-Water-Carbon  Dioxide  Simulation  System 
4.2.1  Calculated  Activity  Coefficients 

The  calculated  activity  coefficients  of  the  chemical 
species  in  the  ascorbic  acid-water-carbon  dioxide  simulation 
system  include  the  activity  coefficients  already  given  in 
Table  4-1  for  the  pure  water-carbon  dioxide  system  and  three 
other  activity  coefficients  associated  with  the  dissociation 
of  ascorbic  acid  given  in  Table  4-5.  Table  4-5  also  gives  the 
corresponding  method  used  to  calculate  these  activity 
coefficients.  In  order  to  calculate  the  variation  of  the 
activity  coefficients  with  the  process  pressure,  the 
compressibility  coefficients  and  the  partial  molar  volumes 
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Table  4-3.  Measured  pH  values  for  pure  water  and  sucrose 
solution 


Process  Pressure, 
P(MPa) 

Average  Measured 
pH  of  Pure  Water  ± 
s 

Average  Measured 
pH  of  20%  Aqueous 
Solution  of 
Sucrose  ±  s 

0 

5.68  ±  0.04 

5.43  ±  0.07 

0.690 

3.55  ±  0.01 

3.51  ±  0.03 

1.379 

3.42  ±  0.00 

3.30  ±  0.01 

2.758 

3.30  ±  0.01 

3.23  ±  0.02 

4.137 

3.23  ±  0.00 

3.21  ±  0.01 

5.861 

3.19  ±  0.00 

3.18  ±  0.01 

Ha:  ^OX*fi20% 
d=0.085 

sd=0.129 

t^  -1.619  <  t0  005=3.365  (for  df=5)  4  Rejection  of  Hypothesis 
(pH  values  are  not  significantly  different  at  0.5%  level.) 
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Table  4-4.  Measured  pH  values  of  water-carbon  dioxide  system 
at  different  process  temperatures  


Process 
rressure f 
P(MPa) 

Measured  pH 

at    Jz  v,   I  S 

Measured  pH 
at  37  C  I  s 

Measured  pH 
at  42  C  ±  s 

0 

5.68  ±  0.04 

5.70  ±  0.14 

5.80  ±  0.00 

0.690 

3.55  ±  0.01 

3.57  ±  0.01 

3.66  ±  0.01 

1.379 

3.42  ±  0.00 

3.44  ±  0.00 

3.51  ±  0.01 

2.069 

3.35  ±  0.01 

3.36  ±  0.01 

3.42  ±  0.00 

2.758 

3.30  ±  0.01 

3.29  ±  0.01 

3.33  ±  0.01 

4. 137 

3.23  ±  0.00 

3.23  ±  0.01 

3.26  ±  0.00 

5.516 

3.19  ±  0.00 

3.20  ±  0.01 

3.24  ±  0.01 

Ho:    %2C=M37C   and  M32C=M42C 

Ha:  Ihn+IHn  and  tHx+Piac 
d^O.010  ;  d2=0.071 

Sd1=0.0l57   ;  sd2=0.0847 

W"1'681  and  2.231  <  t0  005=  3.  143  (for  df=6)  Rejection  of 
Hypothesis  (pH  values  are  not  significantly  different  at  0.5% 
level . ) 
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have  to  be  known  for  ascorbic  acid  and  its  dissociation  ions. 
Since  these  values  were  not  found  in  the  thermodynamic 
literature  and  only  very  high  process  pressures  can  alter 
their  values  significantly,  the  activity  coefficients  for 
ascorbic  acid  and  its  two  uni-  and  bi-  valent  ions  were 
assumed  to  stay  constant  over  the  experimental  pressure  range. 


Table  4-5.  Calculated  activity  coefficients  for  the  ascorbic 
acid-water-carbon  dioxide  system  


Chemical  Species 

Calculated 
Activity 
Coefficient,  y 

Estimation  Method 

C6H8°6 

0.190 

UNIFAC 

C6H706" 

0.760 

Debye-Huckel 

|  W>6" 

0.360 

Debye-Huckel 

4.2.2  Comparison  of  Calculated  and  Experimental  pH  values 

The  pH  data  along  with  the  pH  values  predicted  by  the 
thermodynamic  model  are  given  in  Figures  4-3  and  4-4.  As 
shown  in  the  figures,  the  predicted  pH  values  are  in  good 
agreement  with  the  measured  pH  values.  Tables  A-2  and  A-3 
compare  the  measured  and  predicted  values  of  pH  for  the 
ternary    ascorbic    acid-water-carbon    dioxide    system.  The 
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concentrations  of  ascorbic  acid  in  the  aqueous  solutions  were 
1.107xl0"3  M  and  5.540xl0"4  M  for  the  data  given  in  Tables  A-2 
and  A-3,  respectively.  Over  the  experimental  pressure  range, 
the  pH  of  ascorbic  acid-water-carbon  dioxide  solution  starts 
at  a  lower  value  because  of  the  presence  of  ascorbic  acid  in 
the  solution,  however,  the  pH  drop  from  0  to  5.516  MPa  was  not 
as  great  as  the  one  observed  in  the  binary  water-carbon 
dioxide  system.  From  the  pH  data  it  may  be  concluded  that  the 
presence  of  ascorbic  acid  had  contributed  to  the  reduced 
solubility  of  carbon  dioxide  in  the  solution. 


4.3  Citric  Acid-Water-Carbon  Dioxide  Simulation  System 
4.3.1  Calculated  activity  coefficients 

The  activity  coefficients  for  the  ternary  system  of 
citric  acid-water-carbon  dioxide  simulation  model  includes  the 
activity  coefficients  already  given  for  the  chemical  species 
in  the  pure  water-carbon  dioxide  system  already  given  and  the 
activity  coefficients  involved  in  the  dissociation  of  citric 
acid.  Table  4-6  gives  the  activity  coefficients  calculated 
for  citric  acid  and  its  three  dissociated  ions  along  with  the 
method  of  estimation. 
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Table  4-6.  Calculated  activity  coefficients  for  the  citric 
acid-water-carbon  dioxide  system  


Chemical  Species 

Calculated 
Activity 
Coefficient,  y 

Estimation  Method 

1  C6H807 

0.262 

UNIFAC 

1  W)/ 

0.760 

Debye-Huckel 

I  C6H607" 

0.360 

Debye-Huckel 

C6H507- 

0.115 

Debye-Huckel 

4.3.2  Comparison  of  calculated  and  experimental  pH  values 

Tables  A-4  and  A-5  compare  the  pH  values  measured 
experimentally  with  the  values  predicted  by  the  thermodynamic 
model  for  different  concentrations  of  citric  acids  for  the 
ternary  simulation  model.  These  pH  values  are  plotted  in 
Figures  4-5  and  4-6  where  the  average  measured  pH  is  indicated 
by  data  points  and  the  pH  predicted  by  the  model  is  shown  as 
solid  lines.  In  general,  the  thermodynamic  model  closely 
predicted  the  actual  pH  of  the  simulation  solution  at 
different  pressures  for  various  concentrations  of  citric  acid. 
As  in  the  case  with  the  ascorbic  acid-water-carbon  dioxide 
system,  the  pH  drop  over  the  experimental  pressure  range  was 
not  as  large  as  the  one  for  the  pure  water-carbon  dioxide 
system.     The  same  argument  can  be  made  that  the  presence  of 
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ascorbic  and/or  citric  acid  may  have  decreased  the  solubility 
of  carbon  dioxide  in  the  solution  and  consequently  the  smaller 
drop  in  the  pH. 

4.4  Ascorbic  Acid-Citric  Acid-Water-Carbon  Dioxide  Simulation 

System 

4.4.1  Calculated  Activity  Coefficients 

The  activity  coefficients  used  to  develop  the  pH  model 
for  the  quaternary  simulation  system  of  ascorbic  acid-citric 
acid-water-carbon  dioxide  were  the  ones  calculated  for  the 
three  previous  simulation  models  given  in  Tables  4-1,  4-5  and 
4-6. 

4.4.2  Comparison  of  Calculated  and  Experimental  pH  Values 
Tables  A-6  through  A-10  give  the  measured  and  predicted 

pH  values  for  the  ascorbic  acid-citric  acid-water-carbon 
dioxide  simulation  system  at  different  concentrations  of 
ascorbic  and  citric  acids.  The  plot  of  experimental  and 
calculated  pH  values  as  functions  of  process  pressure  are 
given  in  Figures  4-7  to  4-10.  As  shown  in  these  figures,  the 
drop  in  pH  over  the  experimental  pressure  range  was  even  less 
than  for  the  three  previous  simulation  models.  In  the 
ascorbic  acid-citric  acid-water-carbon  dioxide  system,  the 
presence  of  both  ascorbic  and  citric  acids  seems  to  limit 
carbon  dioxide  solubility  even  more  than  in  the  other 
simulation  systems. 
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4.5  Single  Strength  Orange  Juice 

The  measured  pH  values  at  different  pressures  of  single 
strength  orange  juice  are  compared  with  the  calculated  pH  of 
the    ascorbic    acid-citric    acid-water-carbon    dioxide  model 
solution  in  Figure  4-11  and  Table  A-ll.      The  pH  of  orange 
juice   remained   essentially   constant   over  the  experimental 
pressure  range,   while  the  ascorbic  acid-citric  acid-water- 
carbon  dioxide  model  showed  a  modest  initial  decrease  in  pH 
with  increasing  pressure.     This  may  be  due  to  the  fact  that 
orange  juice  has  a  complex  composition,   and  the  presence  of 
buffer  systems  in  the  juice  can  limit  the  solubility  of  carbon 
dioxide  in  the  juice  and  maintain  the  concentration  of  the 
hydrogen  ions  at  a  fixed  level.  The  presence  of  other  organic 
acids    in    citrus    juices    has    been    indicated    by  several 
investigators.  Nelson  (1928)  reported  traces  of  malic  acid  in 
orange  juice,  and  Hartman  and  Hillig  (1934)  found  appreciable 
amounts  of  it  in  Florida  Valencia  oranges.  Investigating  the 
organic  acids  in  grapefruit,  Menchikowsky  and  Popper  (1932) 
found  that  although  the  organic  acids  of  the  juice  consist 
largely   of   citric   acid    (98.72%),    they   contain   also  malic 
(1.00%),     oxalic      (0.23%),      and     tartaric      (0.05%)  acids 
(percentage  of  the  total  amount  of  acid) . 

4.6  Model  Sensitivity  Analysis 

Tables  4-7  and  4-8  show  the  sensitivity  of  the  pH  models 
to  the  10%  variation  (decrease  or  increase)  in  the  individual 
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activity  coefficients  and  the  reaction  equilibrium  constants. 
The  activity  coefficients  of  hydrogen  ions,  ascorbic  acid  and 
its  univalent  ion,  citric  acid  and  its  uni-  and  bi-valent  ions 
had  a  greater  effect  on  the  predicted  pH  value.  Among  the 
equilibrium  constants,  only  the  first  dissociation  constant 
of  ascorbic  acid  (K4)  and  the  first  dissociation  constant  of 
citric  acid  (K6)  were  found  to  have  some  effect  on  the 
predicted  pH  value.  Increasing  the  values  for  all  the  activity 
coefficients  and  equilibrium  constants  in  the  model  at  the 
same  time  resulted  in  a  3.5%  increase  in  the  predicted  pH 
value.  However,  simultaneous  increase  of  the  activity 
coefficients  and  equilibrium  constants  may  not  indicate  any 
important  conclusions  since  there  are  both  anions  and  cations 
present  in  the  simulation  systems. 

The  major  variables  which  showed  a  significant  effect  on 
the  value  of  predicted  pH  in  the  simulation  models  were  the 
solubility  of  carbon  dioxide  in  water  in  the  water-carbon 
dioxide  system,  and  the  concentrations  of  ascorbic  and  citric 
acids  in  the  three  other  simulation  systems.  In  the  water- 
carbon  dioxide  solution,  doubling  the  solubility  of  carbon 
dioxide  in  the  model  resulted  in  a  3.2  and  4.7%  decrease  in 
the  predicted  pH  value  at  0  and  5.516  MPa  process  pressures, 
respectively.  In  the  quaternary  model,  at  atmospheric 
pressure,  doubling  the  concentrations  of  ascorbic  and  citric 
acids  resulted  in  a  0.7  and  3.9%  drop  in  the  predicted  pH.  The 
concentration  of  citric  acid  was  found  to  have  a  much  greater 
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effect  on  the  predicted  pH  than  the  concentration  of  ascorbic 
acid  for  the  two  ternary  as  well  as  the  quaternary  models, 
which  was  supported  by  the  experimental  data. 

4.7  Discussion  of  Results 
The  mathematical  models  developed  for  the  four  simulation 
systems  accurately  predicted  the  pH  variation  of  these 
solutions  as  functions  of  process  pressure  and  the 
concentrations  of  ascorbic  and  citric  acids.  From  the 
expressions  developed  for  the  variation  of  the  equilibrium 
constants  and  the  activity  coefficients  as  a  function  of 
process  pressure,  it  was  concluded  that  pressure  has  a 
significant  effect  only  at  very  high  values,  and  that  it  is 
the  solubility  of  carbon  dioxide  in  water  which  has  a  major 
effect  on  the  pH  of  the  solution.  For  the  sugar  concentration 
and  process  temperatures  tested,  it  was  shown  that  the  pH  was 
not  a  strong  function  of  sugar  concentration  or  process 
temperature . 

The  predicted  and  measured  pH  values  for  the  binary 
water-carbon  dioxide  system  showed  the  greatest  drop  in 
solution  pH  from  5.70  to  3.20  over  the  pressure  range  of  0  to 
5.52  KPa  at  a  temperature  of  37°C.  This  implicated  that 
substantially  greater  quantities  of  carbon  dioxide  were 
dissolved  to  produce  carbonic  acid  in  the  binary  water-carbon 
dioxide  system  than  in  the  other  simulation  systems  resulting 
in    the    lowering    of    pH   values    in    the    system    during  the 
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pressurization  process.  The  pH  drop  from  the  initial  value  at 
atmospheric  pressure  in  the  three  simulation  systems  of 
ascorbic  acid-water-carbon  dioxide,  citric  acid-water-carbon 
dioxide,  and  ascorbic  acid-citric  acid-water-carbon  dioxide 
was  found  to  be  considerably  less  than  in  the  binary  system 
mainly  because  of  the  presence  of  ascorbic  acid  and  citric 
acid  in  the  solution  initially,  which  presumably  limited  the 
solubility  of  carbon  dioxide  and  the  subsequent  production  of 
carbonic  acid. 

This  behavior  of  the  system  pH  was  even  more  pronounced 
in  the  experiments  performed  with  single  strength  orange  juice 
over  the  pressure  range  of  0  to  5.57  MPa.  The  pH  of  the 
orange  juice  dropped  by  only  3.2%  whereas  for  the  pure  water- 
carbon  dioxide  solution  the  pH  drop  was  43.86%  over  the  same 
pressure  range.  From  Figure  4-11  which  shows  both  the 
measured  pH  of  the  orange  juice  and  the  pH  values  predicted  by 
the  ascorbic  acid-citric  acid-water-carbon  dioxide  model,  it 
may  be  concluded  that  the  complex  composition  of  orange  juice, 
or  citrus  juice  in  general,  and  the  existing  buffer  systems 
are  a  major  factor  in  the  low  sensitivity  of  pH  with  process 
pressure.  The  small  pH  drop  also  indicates  that  cloud 
stability  and  pectinesterase  inactivation  during  the  high 
pressure  carbon  dioxide  process  may  not  be  entirely  due  to  low 
pH  inactivation  of  the  enzyme,  and  that  other  mechanisms  of 
inactivation  may  be  responsible  for  cloud  stability. 
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Table  4-7.  Sensitivity  analysis  of  the  pH  model  at  atmospheric 
pressure  and  37°C  with  the  variation  in  activity  coefficients 


Activity 
coefficient 

-10% 
Variation 

0  2:  Va  r  i  a    "i  on 

Variation 

H+ 

3.365 

3.343 

3.324 

OH" 

3.343 

3.343 

3.343 

co2 

3.343 

3.343 

3.343 

H2C03 

3.344 

3.343 

3.343 

HC03" 

3.343 

3.343 

3.343 

co32" 

3.343 

3.343 

3.343 

C6H806 

3.346 

3.343 

3.341 

C6H706* 

3.341 

3.343 

3.345 

1  W^" 

3.343 

3.343 

3.343 

1  C6H80? 

3.362 

3.343 

3.326 

C6H70/ 

3.327 

3.343 

3.358 

3.341 

3.343 

3.346 

C6H5073- 

3.343 

3.343 

3.343 
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Table  4-8.  Sensitivity  analysis  of  the  pH  model  at  atmospheric 
pressure  and  37°C  with  the  variation  in  the  equilibrium 
constants 


Equilibrium 
constant 

-10% 
Variation 

0%  Variation 

+  10% 
Variation 

K 

3.343 

3.343 

3.343 

Ki 

3.343 

3.343 

3.343 

1  K2 

3.344 

3.343 

3.343 

K3 

3.343 

3.343 

3.343 

** 

3  .  346 

3.343 

3.341 

v 

7  TAT 

1     *S  A  *3 

J  •  j4  j 

J  .  J4  J 

!  K6 

3.362 

3.343 

3.326 

l  K? 

3.343 

3.343 

3.343 

1  *8 

3.343 

3.343 

3.343 

CHAPTER  5 

CONCLUSIONS     AND     SUGGESTIONS  FOR  FUTURE  STUDY 
5.1  Conclusions 
Thermodynamic  correlations  were  successfully  applied  to 
a  set  of  four  citrus  simulation  systems  to  predict  their  pH 
variation  when  processed  with  supercritical  carbon  dioxide, 
with  the  following  conclusions  drawn  from  this  study: 

1.  Mathematical  expressions  can  be  developed  to  model 
the  behavior  of  food  systems    such  as  fruit  juices  through  a 
thermodynamic  approach.  All  of  the  pH  models  developed  for 
the  four  simulation  solutions  of  pure  water-carbon  dioxide, 
ascorbic  acid-water-carbon  dioxide,  citric  acid-water-carbon 
dioxide  and  ascorbic  acid-citric  acid-water-carbon  dioxide, 
accurately  predicted  the  pH  values  which  were  measured 
experimentally.  Over  the  experimental  pressure  range,  the 
lowering  of  the  pH  of  the  simulation  systems  was  mostly  due 
to  the  solubility  of  carbon  dioxide  and/or  the 
concentrations  of  acids  in  the  solution  rather  than  the 
variation  of  the  equilibrium  constants  or  the  activity 
coefficients  with  pressure. 

2.  The  pH  values  for  both  the  quaternary  simulation 
system  and  single  strength  orange  juice  were  found  to  be 
weak  functions  of  process  pressure  but  very  strong  functions 
of  system  composition.  The  pH  of  the  orange  juice  was  less 
sensitive  to  process  pressure  than  the  pH  of  the  ascorbic 
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acid-citric  acid-water-carbon  dioxide  system.  This  may  be 
due  to  the  complex  composition  of  citrus  juices  which,  with 
their  natural  buffer  systems,  can  limit  the  solubility  of 
carbon  dioxide    and  consequently  limit  the  production  of 
carbonic  acid  in  the  solution.     This  implication  was  deduced 
from  the  observation  that  the  pH  drop  in  the  binary  system 
of  water-carbon  dioxide  with  a  relatively  simple  composition 
was  43.86%,  whereas  the  pH  drop  for  the  single  strength 
orange  juice  was  only  3.2%  over  the  same  pressure  range. 
Generally  as  the  composition  of  the  simulation  systems 
became  more  complex,  less  carbon  dioxide  was  dissolved  in 
the  solutions.  Since  this  limits  the  generation  of  carbonic 
acid,  the  observed  pH  drop  was  smaller. 

3.  The  pH  values  predicted  by  the  models  for  the  four 
simulation  systems    can  provide  information  about  pH 
variations  of  not  only    citrus  juices  but  also  citrus 
beverages  during  the  carbonation  process.  From  the  pH 
studies  with  the  four  simulation  systems  it  was  concluded 
that  the  pH  drop  for  a  citrus  beverage  during  the 
carbonation  process  can  be  composition  dependent  and/or 
pressure  dependent. 

4.  The  small  pH  drop  in  orange  juice  and  the  ternary 
and  quaternary  models  did  not  support  the  hypothesis  that 
pectinesterase  was  inactivated  by  low  pH  in  the  high 
pressure  carbon  dioxide  process.  Arreola  (1990)  reported 
100%  inactivation  of  pectinesterase  in  orange  juice  treated 
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with  supercritical  carbon  dioxide  at  pressures  between  7.4 
and  34  MPa  and  temperatures  of  35  to  60°C.  However,  the  pH 
behavior  of  orange  juice  was  not  measured  during  those 
experiments  and  was  unknown.     The  mechanism  of 
pectinesterase  inactivation  and  the  subsequent  cloud 
stability  in  the  high  pressure  carbon  dioxide  treatment  of 
citrus  juices  may  only  be  partly  due  to  lowering  of  the  pH 
if  at  all;  it  may  be  strictly  a  pressure-related  phenomena. 

5.2  Suggestions  For  Future  Study 
Future  works  in  the  area  of  citrus  cloud  stability  in 
the  supercritical  carbon  dioxide  process  can  focus  on  the 
mechanism  of  pectinesterase  inactivation  and  its  relation  to 
the  pH.  More  studies  on  the  inactivation  of  pectinesterase 
are  required  to  not  only  find  out  what  the  mechanism  is,  but 
also  to  be  able  to  quantify  the  percentage  of  inactivation 
of  pectinesterase.  The  same  way    thermodynamics  and  reaction 
equilibria  were  applied  in  this  study  to  find  mathematical 
expressions  for  the  pH  variation  of  the  simulation  systems, 
principles  of    reaction  equilibria  and  enzyme  kinetics  can 
be  used  together  to  quantitatively  model  cloud  stability  in 
citrus  juices  as  a  function  of  the  supercritical  process 
conditions. 


APPENDIX  A 

TABLES  OF  PREDICTED  AND  MEASURED  PH  VALUES 
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Table  A-l.  Predicted  and  measured  pH  values  at  different 
process  pressures  for  the  pure  water-carbon  dioxide 
simulation  system  


Process  Pressure, 
P(MPa) 

Average 
Experimental  pH  ± 
s 

Predicted  pH 

0 

5.75  ±  0.05 

5.59 

0.3448 

4.59  ±  0.09 

4.79 

0. 6895 

3.62  ±  0.10 

3.90 

1.  379 

3.48  ±  0.06 

3.41 

2.758 

3.32  ±  0.04 

3.28 

4.138 

3.27  ±  0.02 

3.21 

5.516 

3.22  ±  0.03 

3.17 

6.  895 

3.17  ±  0.04 

3.14 

10.343 

3 . 15  ±  0. 04 

3.12 

13  . 790 

3.15  ±  0.03 

3.11 

17 .238 

3.14  ±  0.04 

3.10 

20. 685 

3.13  ±  0.04 

3.09 

24 . 133 

3.13   +   0  04 

27.580 

3.12  ±  0.03 

3.08 

31.028 

3.10  ±  0.04 

3.07 

34.475 

3.08  ±  0.04 

3.06 

Table  A-2 .  Comparison  of  measured  and  predicted  pH  values  of 
the  ascorbic  acid-water-carbon  dioxide  system  containing 
1.107x10    M  ascorbic  acid 


Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.74  ±  0.15 

3.85 

0.690 

3.47  ±  0.08 

3.51 

1.379 

3.27  ±  0.04 

3.39 

2.758 

3.25  ±  0.05 

3.28 

4. 137 

3.20  ±  0.01 

3.21 

5.516 

3.13  ±  0.03 

3.16 

Table  A-3 .  Comparison  of  measured  and  predicted  pH  values  of 
the  ascorbic  acid-water-carbon  dioxide  system  containing 
5.540x10     M  ascorbic  acid 

Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.90  ±  0.08 

3.98 

0.690 

3.45  ±  0.06 

3.53 

1.379 

3.35  ±  0.11 

3.40 

2.758 

3.26  ±  0.03 

3.26 

4.137 

3.19  ±  0.01 

3.19 

5.516 

3.15  ±  0.03 

3.15 
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Table  A-4.  Comparison  of  measured  and  predicted  pH  for  the 
citric  acid-water-carbon  dioxide  system  with  6.08x10"  M 
citric  acid 




Process  Pressure, 

p (MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.32  ±  0.04 

3.35 

0.414 

3.23  ±  0.00 

3.32 

0.690 

3.20  ±  0.05 

3.28 

1.  379 

3.16  ±  0.07 

3.22 

2.758 

3.12  ±  0.02 

3.14 

4.137 

3.08  ±  0.06 

3.10 

5.516 

3.05  ±  0.03 

3.06 

Table  A-5.  Comparison  of  measured  and  predicted  pH  for  the 
citric  acid-water-carbon  dioxide  system  with  3.04xl0"4  M 
citric  acid 

Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.52  ±  0.04 

3.48 

0.690 

3.31  ±  0.07 

3.36 

1.379 

3.26  ±  0.01 

3.28 

2.758 

3.15  ±  0.03 

3.18 

4.137 

3.12  ±  0.04 

3.12 

5.516 

3.10  ±  0.00 

3.09 
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Table  A-6.  Comparison  of  measured  and  predicted  pH  for  the 
ascorbic^  acid-citric  acid-water-carbon  dioxide  system  with 
5.54x10    M  ascorbic  acid  and  6.08x10    M  citric  acid. 


Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

■i 

Predicted  pH 

0 

3.33  ±  0.09 

3.35 

0.276 

3.29  ±  0.05 

3.32 

0.690 

3.25  ±  0.01 

3.28 

1 .  379 

3.20  ±  0.04 

3  .22 

2.758 

3.12  ±  0.06 

3.14 

4.137 

3.09  ±  0.05 

3.10 

jl  5.516 

3.07  ±  0.02 

3.06 

Table  A-7.  Comparison  of  measured  and  predicted  pH  for  the 
ascorbic  3acid-citric  acid-water-carbon  dioxide  system  with 
1.107x10     M  ascorbic  acid  and  6.08x10"   M  citric  acid 

Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.32  ±  0.06 

3.34 

0.690 

3.30  ±  0.02 

3.27 

1.379 

3.17  ±  0.05 

3.21 

2.758 

3.12  ±  0.08 

3.13 

4.137 

3.11  ±  0.01 

3.09 

5.516 

3.06  ±  0.03 

3.06 

113 


Table  A-8.  Comparison  of  measured  and  predicted  pH  for 
ascorbic  3acid-citric  acid-water-carbon  dioxide  system  with 
1.107x10"   M  ascorbic  acid  and  3.04x10    M  citric  acid 


Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.44  ±  0.03 

3.46 

0. 690 

3.27  ±  0.05 

3.35 

1 . 379 

3.20  ±  0.11 

3.26 

2.758 

3.13  ±  0.08 

3.17 

4.137 

3.09  ±  0.04 

3.12 

5.516 

3.10  ±  0.07 

3.08 

Table  A-9.  Comparison  of  measured  and  predicted  pH  values 
for  the  ascorbic  acid-citric  acid-water-carbon  dioxide 
system  with  5.54x10"   M  ascorbic  acid  and  3.04xl0*4  M  citric 
acid. 

Process  Pressure, 
P(MPa) 

Measured  pH  ±  s 

Predicted  pH 

0 

3.46  ±  0.13 

3.47 

0.345 

3.31  ±  0.05 

3.41 

0.690 

3.28  ±  0.02 

3.36 

1.379 

3.22  ±  0.04 

3.27 

2.758 

3.14  ±  0.02 

3.18 

4.137 

3.12  ±  0.05 

3.13 

5.516 

3.08  ±  0.01 

3.09 
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Table  A-10.  Comparison  of  measured  and  predicted  pH  of 
single  strength  orange  juice  


Process  Pressure . 
P(MPa) 

0 

3.75  ±  0.13 

3.80 

0.690 

3.73  ±  0.05 

3.47 

1.379 

3.68  ±  0.06 

3.34 

2.758 

3.66  ±  0.02 

3.23 

4.137 

3.65  ±  0.03 

3.16 

5.516 

3.63  ±  0.00 

3.11 

APPENDIX  B 


FORTRAN  PROGRAM  TO  SOLVE  FOR 
EQUILIBRIUM  ACTIVITIES  USING 
NEWTON-RAPHSON  ITERATION  METHOD 
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******************************************************* 

*  MODEL  FOR  THE  PREDICTION  OF  PH  AS  A  FUNCTION  OF  PRESSURE  AND 

*  TEMPERATURE  FOR  THE  SIMULATION  SYSTEMS  OF  A) PURE  WATER-CARBON 

*  DIOXIDE   ;   B) ASCORBIC  ACID-CITRIC  ACID-WATER-CARBON  DIOXIDE 
**************************************************************** 

*  THIS  PROGRAM  CALCULATES  PH  OF  THE  SIMULATION  SYSTEM  FOR  A  SPE- 
CIFIED PROCESS  TEMPERATURE  AND  PRESSURE  ,  AND  CONCENTRATIONS 
OF  ASCORBIC  ACID  AND  CITRIC  ACID  USING  THE  NEWTON-RAPHSON 
ITERATION  METHOD  BY  SOLVING  A  SET  OF  SIMULTANEOUS  EQUATIONS  FOR 
THE  EQUILIBRIUM  CONCENTRATIONS  OF  THE  CHEMICAL  SPECIES  PRESENT 
IN  THE  SYSTEM.  THE  PH  OF  THE  SOLUTION  IS  CALCULATED  BY  TAKING 
THE  NEGATIVE  LOG  OF  THE  CONCENTRATION  OF  HYDROGEN  IONS  TO 
WHICH  THE  ITERATION  VARIABLES  CONVERGE. 

**************************************************************** 

*  DACTI=ACTIVITY  COEFFICIENT  OF  SPECIES  I 

*  DX (I) =EQUI LIBRIUM  CONCENTRATION  OF  SPECIES  I  (MOLES/LITER) 


* 

SPECIES 

1= 

H+ 

IONS 

* 

SPECIES 

2= 

OH- 

IONS 

* 

SPECIES 

3= 

C02 

(CARBON  DIOXIDE) 

* 

SPECIES 

4= 

H2C03 

(CARBONIC  ACID) 

* 

SPECIES 

5= 

HC03- 

IONS 

* 

SPECIES 

6= 

C03  — 

IONS 

* 

SPECIES 

7= 

C6H706- 

IONS 

* 

SPECIES 

8= 

C6H606 — 

IONS 

* 

SPECIES 

9= 

C6H707- 

IONS 

* 

SPECIES 

10= 

C6H607 — 

IONS 

* 

SPECIES 

11= 

C6H507  

IONS 

**************************************************************** 
*  DKW=DISSOCIATION  CONSTANT  FOR  DISSOCIATION  OF  WATER 
DKJ=DISSOCIATION  CONSTANT  FOR  THE  REACTION  J 


REACTION  2= 
REACTION  3= 
REACTION  4= 


REACTION  5= 
REACTION  6= 
REACTION  7= 
REACTION  8= 
REACTION  9= 


DISSOCIATION  OF  CARBONATE  ION 

DISSOCIATION  OF  CARBONIC  ACID  TO  CARBONATE  ION 
DISSOCIATION  OF  CARBONIC  ACID  TO  CARBON  DIOXIDE 
AND  WATER 

DISSOCIATION  OF  ASCORBIC  ACID  TO  C6H706-  ION 

DISSOCIATION  OF  C6H706-  TO  C6H606 — 

DISSOCIATION  OF  CITRIC  ACID  TO  C6H707- 

DISSOCIATION  OF  C6H707-  TO  C6H607 — 

DISSOCIATION  OF  C6H607 —  TO  C6H507  

********************************************** ****************** 

*  DSOL=SOLUBILITY  OF  CARBON  DIOXIDE  IN  WATER 

*  DFUN2=INTERPOLATION  FUNCTION  FOR  CALCULATING  THE  SOLUBILITY 
**************************************************************** 

INTEGER  SYS,I,J,L,K,Q 

DOUBLE  PRECISION  DACT1 , DACT2 , DACT3 , DACT4 , DACT5 , DACT6 
DOUBLE  PRECISION  DKW , DK2 , DK3 , DK4 , DK5 , DH , DFH , DDFH , DNEWH 
DOUBLE  PRECISION  DXH20, DA, DB, DC, DD, DE, DF, DP, DPH, DS , DETM 
DOUBLE  PRECISION  DX ( 11) , B ( 11) , C ( 11 , 11) , Bl ( 11) , DXX ( 11) , DK6 
DOUBLE  PRECISION  DTNRM, DXCIT , DXASC , DACT7 , DACT8 , DACT9 
DOUBLE  PRECISION  DACT10 , DACT11 , DK7 , DK8 , DK9 , DK10 , DK11 
DOUBLE  PRECISION  PD,CC 
DOUBLE  PRECISION  P(90) , T (90) , S (90) 
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DOUBLE  PRECISION  DEN , DSOL , DTEMP , DFUN2 

PRINT* , • ENTER  SIMULATION     MODEL  SYSTEM:' 

PRINT* ,'1     FOR  PURE  WATER-CARBON  DIOXIDE  SYSTEM' 

PRINT*,' OR  2  FOR  ASCORBIC-CITRIC-WATER-CARBON  DIOXIDE 

SYSTEM ' 
READ* , SYS 

**************************************************** 

*  INTERPOLATION  ROUTINE  FOR  THE  CALCULATION  OF  CARBON  DIOXIDE 

*  SOLUBILITY  AT  THE  SPECIFIC  TEMPERATURE  AND  PRESSURE 
*********************************************** 

************************************************************** 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


P(N),N=1,9)  /9*25.0D+00/ 
P(N) ,N=10, 18)  /9*50.0D+00/ 
P(N) ,N=19,27)  /9*75.0D+00/ 
P(N) ,N=28,36)  /9*100.0D+00/ 
P(N) ,N=37,45)  /9*150.0D+00/ 
P(N) ,N=46,54)  /9*200.0D+00/ 
P(N) ,N=55,63)  /9*300.0D+00/ 
P(N) ,N=64, 72)  /9*400.0D+00/ 
P(N) ,N=73,81)  /9*500.0D+00/ 
P(N) ,N=82,90)  /9*700.0D+00/ 


************************************************************* 


T(N) ,N=1,9)   /12, 18, 25, 31. 04, 35, 40 ,50, 75, 100/ 

T(N) ,N=10,18)  /12, 18, 25, 31. 04, 35, 40, 50, 75, 100/ 

T(N) ,N=19,27)  /12, 18, 25, 3 1.04, 35, 40, 50, 75, 100/ 

T(N) ,N=28,36)  /12 , 18 , 25 , 31 . 04 , 35 , 40 , 50 , 75 , 100/ 

T(N) ,N=37,45)  /12 , 18 , 25 , 31 . 04 , 35 , 40 , 50 , 75 , 100/ 

T(N) ,N=46,54)  /12 , 18 , 25 , 31 . 04 , 35 , 40 , 50 , 75 , 100/ 

T(N) ,N=55,63)  /12 , 18 , 25 , 3 1 . 04 , 35 , 40 , 50 , 75 , 100/ 

T(N) ,N=64,72)  /12 , 18 , 25 , 31 . 04 , 35 , 40 , 50 , 75 , 100/ 

T(N) ,N=73,81)  /12, 18, 25, 31. 04, 35, 40, 50, 75, 100/ 

T(N) ,N=82,90)  /12, 18, 25, 31. 04, 35, 40, 50, 75, 100/ 


********************************************** ******************* 
DATA(S(N) , N=l, 9) /6. 88, 3. 86, 4. 59, 2. 80, 2. 56, 2. 30, 1.92, 1.35, 1.06/ 
DATA (S (N) ,N=10, 18) /7. 03, 6. 33, 5. 38, 4. 77, 4. 39, 4. 02, 3. 41, 2. 49, 2. 01/ 
DATA (S (N) ,N=19,27)/7.18,6.69,6.17,5.8,5.51,5.1,4.45,3.37,2.82/ 
DATA (S (N) ,N=28,36)/7.27,6.72,6.28,5.97,5.76,5.5,5.07,4.07,3.49/ 
DATA(S(N) ,N=37, 45) /7. 59, 7. 07, 6. 60, 6. 25, 6. 03, 5. 81, 5. 47, 4. 86, 4. 49/ 
DATA(S(N) ,N=46,54)/7.6,7.21,6.91,6.48,6.29,6.28,5.76,5.27,5.08/ 
DATA (S (N) ,N=55,63)/7.86,7.35,7.23,6.88,6.68,6.59,6.2,5.83,5.84/ 
DATA (S (N) ,N=64,72)/8.12,7.77,7.54,7.27,7.06,6.89,6.58,6.3,6.4/ 
DATA(S(N) ,N=73,81)/8.38,8.19,7.85,7.65,7.51,7.26,6.96,6.4,6.96/ 
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DATA(S(N) , N=82, 90) /8. 9, 8. 61, 8. 47, 8. 42, 8. 34, 7. 93, 7. 58, 7. 43, 7. 61/ 

**************************************************** 

***************************************************************** 

IF   (SYS.EQ.l)  THEN 

PRINT* ,  •  • 

101     PRINT*, 'ENTER  PRESSURE  (PSI)=' 
READ* , DP 

IF   (DP.LT.0.0D+00)  THEN 

STOP 

ELSE 

GOTO  110 
ENDIF 

110       PRINT* ,• ENTER  TEMPERATURE  (C)=' 
READ* , DTEMP 

PRINT* ,  1  ' 

DP=DP/14.7D+00 
DO  10  Q=l,l 

N=90 

M=9 

DX(l)=1.0D-4 
DO  780  1=2,6 

DX(I)=1.0D-4 
780  CONTINUE 

IF   (DP.LE.24.0D+00)  THEN 

DSOL=( (0.5942D+00) * ( DP-1 . 0) +0 . 696D+00) /2.4D+01 
ELSE 

DSOL=DFUN2 (DP, DTEMP, N,M, P, T, S) *2 . 273D-01 
ENDIF 

********************************************************* 
DKW=(10.0D+00) **( (-13.677D+00)+(0.00036D+00) *DP) 
DK2= (10. OD+00) **( (0.0005D+00) *DP-10 . 2515D+00) 
DK3=( 10. OD+00) **( (0.0005D+00) *DP-6 . 3662D+00) 
DK4=(10. OD+00) **( (-6.355D+00)+(0.00048D+00) *DP) 

****************************************************** *********** 

IF(DP.LE.7.0D+00)  THEN 

DACT1= ( 0 . 0005D+00*DEXP (DP-1 . OD+00) ) 

ELSE 

DACTl=0.830D+00 
ENDIF 

DACT2= (10. OD+00) ** ( (-0 . 11926D+00) - (2 . 575D-06) *DP*DP+ ( 1 . 03D-04 ) *DP 
DACT3=1. OD+00 

DACT4=( 10. OD+00) ** ( (+0 . 0005D+00) - ( 0 . 000479D+00) *DP) 
DACT5= (10 . OD+00) ** ( (-0 . 11068D+00) + (1 . 421D-05) *DP- (3 . 55D-07) *DP*DP 


DACT6= (10. OD+00) **( (-0 . 44362D+00) - ( 6 . 661D-06) *DP*DP+ ( 1 . 332D-05) 
+  *DP) 

****'k**************************1c***1ciflc**********1c**1c**1c*********1c 
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DO  20  J=l,100 

DO  500  1=1,6 

DO  600  K=l,6 

C(I,K)=0.0D+00 
600  CONTINUE 
500  CONTINUE 


C(2 

rl] 

1=1. 0D+00 

C(3 

=-DACTl *  DACT2 *  DX ( 2 ) 

C(4 

r  x; 

=-DACTl*DACT6*DX(6) 

C(5 

r-l'l 

=-DACTl*DACT5*DX(5) 

C(l 

-3; 

=1.0D+00 

C(l 

-4] 

=1.0D+00 

C(l 

=1.0D+00 

C(l, 

6; 

=1.0D+00 

C(2, 

-2; 

=-1. 0D+00 

C(2; 

5; 

=-1. 0D+00 

C(2, 

6, 

=-2 . 0D+00 

C(3, 

2; 

=-DACTl*DACT2*DX(l) 

C(4, 

5; 

=DACT5*DK2 

C(4, 

6; 

=-DACTl*DACT6*DX ( 1) 

C(5, 

4; 

=DACT4*DK3 

C(5( 

5; 

=-DACTl*DACT5*DX(l) 

C(6, 

3; 

=-DACT3 

C(6( 

4; 

=DACT4*DK4 

****************************************** 

B(l)=DX(3)+DX(4)+DX(5)+DX(6)-DSOL 
B(2)=DX(1)-DX(2)-DX(5)-2.0D+00*DX(6) 
B(3)=-DACT1*DACT2*DX(1) *DX(2)   +  DKW 
B(4)=-DACT1*DACT6*DX(1) *DX(6)   +  DACT5*DK2*DX (5) 
B(5)=-DACT1*DACT5*DX(1) *DX(5)   +  DACT4*DK3*DX (4 ) 
B(6)=-DACT3*DX(3)   +  DACT4*DK4*DX(4) 
***************************************************************** 

N=6 

*********************************************** 

CALL  INVDET (C, N, DTNRM , DETM) 
DO  700  K=l,6 

Bl(K)=0.0D+00 
DO  710  L=l, 6 

B1(K)=C(K,L)*B(L)+B1(K) 
710  CONTINUE 
700  CONTINUE 
*********************************************** 
DO  720  L=l, 6 

DXX(L)=DX(L) -B1(L) 
720  CONTINUE 

IF(DXX(2) .GT.0D+00.AND.DXX(3) . GT. 0 . OD+00 . AND. DXX (4 ) .GT. 
+       0. OD+00. AND. DXX (5) .GT. 0. OD+00. AND. DXX(6) . GT. 0 . OD+00 . AND. 

+       DXX(l) .GT.O. OD+00. AND. DABS (DXX(l)-DX(l) ) .LE.1.0D-9)  THEN 
IF(DXX(1) .GT.O. OD+00)  THEN 
PH=-DLOG10(DXX(1) *DACT1) 
PRINT*  ,  •  PH=  '  ,  PH 

PRINT* ,  1  1 


ENDIF 
GOTO  101 
ENDIF 

DO  22  1=1,6 

DX(I)=DXX(I) 
22  CONTINUE 
20  CONTINUE 

********************************************************** 
10  CONTINUE 

************************************************************** 

ELSEIF   (SYS.EQ.2)  THEN 
************************************************************** 
************************************************************** 

DO  11  Q=l,100 

PRINT* ,  1  • 

PRINT* ENTER  PRESSURE  (PSI)=« 
READ* , DP 

PRINT* ,  1  • 

IF(DP.EQ.-1.0D+00)   GOTO  5000 

DXCIT=3.04D-04 

DXASC=1.107D-03 

DTEMP=37.0D+00 

DP=DP/14.7D+00 

N=90 

M=9 

DX(l)=1.0D-04 

DX(2)=1.0D-04 

DX(3)=0.1D-04 

DX(4)=0.1D-04 

DX(5)=0.1D-04 

DX(6)=0.1D-04 

DX(7)=4.0D-04 

DX(8)=4.0D-04 

DX(9)=4.0D-04 

DX(10)=4.0D-04 

DX(ll)=4.0D-04 

IF   (DP.LE.24.0D+00)  THEN 

DSOL=( (0.5942D+00) * (DP-1 . 0) +0 . 696D+00) /2.4D+01 
ELSE 

DSOL=DFUN2 (DP, DTEMP, N,M, P, T, S) *2.273D-01 
ENDIF 

***************************************^AAJk^^AAAitA#jfcic#jlcAA^^^A^AJkA 

DKW= (10. OD+00) ** ( (-13 . 677D+00) + (0. 00036D+00) *DP) 

DK2=(10.0D+00) **( (0.0005D+00) *DP-10 . 2515D+00) 

DK3=(10. OD+00) **( (0.0005D+00) *DP-6 . 3662D+00) 

DK4=(10. OD+00) **( (-6.355D+00)+(0.00048D+00) *DP) 

DK5=7.94D-05 

DK6=1.62D-12 

DK10=1.73D-05 

DKll=4.02D-07 

A************************************^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
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DACTl=0.830D+00 

DACT2=(10.0D+00) ** ( (-0 . 11926D+00) - (2 . 575D-06) *DP*DP+ ( 1 . 03D-04 ) *DP 

) 

DACT3=1.0D+00 

DACT4=(10.0D+00) ** ( (+0 . 0005D+00) - (0 . 000479D+00) *DP) 

DACT5=(10.0D+00)**( (-0 . 11068D+00) + ( 1 . 421D-05) *DP- (3 . 55D-07) *DP*DP 
) 

DACT6=(10.0D+00) ** ( (-0 . 44362D+00) - (6 . 661D-06) *DP*DP+(1. 332D-05) 
+  *DP) 

DACT7=0.76D+00 

DACT8=0.36D+00 

DACT9=0.76D+00 

DACT10=0.36D+00 

DACTll=0.115D+00 
***************************************************** 

DO  21  J=l,200 

DO  501  1=1,11 

DO  601  K-1,11 

C(I,K)=0.0D+00 
601  CONTINUE 
501  CONTINUE 


c 

(2 

pl 

|=1.0D+00 

c 

(3 

,1, 

|=-DACT1*DACT2*DX(2) 

c 

(4 

f  1] 

|=-DACT1*DACT6*DX(6) 

c 

(5 

rl] 

|=-DACT1*DACT5*DX(5) 

c 

(1 

r3] 

|=1.0D+00 

c 

:i 

r4] 

1=1. OD+00 

c 

;i 

=1.0D+00 

c 

;i 

-6; 

=1.0D+00 

c 

[2 

r2] 

=-1. 0D+00 

c 

(2 

r5; 

=-1. OD+00 

c 

[2 

=-2 . OD+00 

c 

!3 

-2; 

=-DACT 1 * DACT2 *  DX ( 1 ) 

c 

;4 

=DACT5*DK2 

c 

;4 

=-DACTl*DACT6*DX(l) 

c 

4; 

=DACT4*DK3 

c 

5; 

=-DACTl*DACT5*DX(l) 

c 

[6, 

3; 

=-DACT3 

c 

[6, 

4; 

=DACT4*DK4 

c 

'J, 

1; 

=DX(7) *DACT1*DACT7 

c 

'J , 

7; 

=DX(1) *DACT1*DACT7 

C| 

'8, 

1; 

=DX(8) *DACT1*DACT8 

C| 

'8, 

7; 

=-DK6*DACT7 

C| 

'8, 

8] 

=DX(1) *DACT1*DACT8 

C( 

2, 

7] 

=-1. OD+00 

C( 

2, 

8) 

=-2. OD+00 

c< 

2, 

9) 

=-1. OD+00 

C( 

2, 

10)=-2.0D+00 

C( 

2, 

ll)=-3.0D+00 
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C ( 9 , 1 ) =-DACT9 *DACT1*DX ( 9 ) 
C(9,9)=-DACT1*DX(1) *DACT9 
C ( 10 , 1) =-DACT10*DACTl*DX (10) 
C ( 10 , 9 ) =DACT9*DK10 
C(10,10)=-DACT10*DACT1*DX(1) 
C(11,1)=-DACT11*DACT1*DX(11) 
C(ll, 10)=DACT10*DK11 
C ( 11 , 11) =-DACTll*DACTl*DX (1) 
***************************************************** 

B(l)=DX(3)+DX(4)+DX(5)+DX(6)-DSOL 

B(2)=DX(l)-DX(2)-DX(5)-(2.0D+00) *DX ( 6) -DX (7 ) - (2 . OD+00) *DX(8) 
+  -DX(9)-2.0D+00*DX(10)-3.0D+00*DX(11) 
B(3)=-DACT1*DACT2*DX(1) *DX(2)   +  DKW 
B(4)=-DACT1*DACT6*DX(1) *DX(6)   +  DACT5*DK2*DX (5) 
B(5)=-DACT1*DACT5*DX(1) *DX(5)   +  DACT4*DK3*DX(4) 
B(6)=-DACT3*DX(3)   +  DACT4*DK4*DX(4) 

B(7)=DACT1*DX(1) *DX(7) *DACT7- (7 . 94D-05) *DXASC*0 . 19D+00 
B(8)=DACT1*DX(1) *DX(8) *DACT8-1 . 62D-12*DX (7) *DACT7 
B(9)=(0.262D+00*7.45D-04) *DXCIT- (DACT9*DACT1*DX ( 1) *DX(9) ) 
B(10)=DACT9*1.73D-05*DX(9)-DACT10*DACT1*DX(1) *DX(10) 
B(11)=DACT10*4.02D-07*DX(10)-DACT11*DACT1*DX(1) *DX(11) 
******************************************************* 
N=ll 

*********************************************** 

CALL  INVDET (C, N, DTNRM, DETM) 
DO  701  K=l,ll 

Bl(K)=0.0D+00 
DO  711  L-1,11 

B1(K)=C(K,L) *B(L)+B1(K) 
711  CONTINUE 
701  CONTINUE 
*********************************************** 

DO  721  L=l,ll 

DXX(L)=DX(L)-B1(L) 
721  CONTINUE 

IF(DXX(2) .GT.0D+00.AND.DXX(3) . GT. 0 . OD+00 . AND. DXX (4 ) .GT. 
+       0. OD+00. AND. DXX (5) . GT. 0 . OD+00 . AND. DXX ( 6) . GT. 0 . OD+00 .AND. 
+       DXX (1) .GT. 0 . OD+00 .AND. DABS (DXX (1) -DX ( 1) ) . LE. 1 . OD-9)  THEN 

IF(DXX(1) .GT.O. OD+00)  THEN 

PH=-DLOG10(DXX(1) *DACT1) 

PRINT* ,  '  PH=  •  ,  PH 

PRINT*,  1  PRESS (PSI)  =  • ,DP*14.7D+00 

PRINT* ,  1  i 

ENDIF 

GOTO  11 

ENDIF 

DO  23  1=1,11 

DX(I)=DXX(I) 
23  CONTINUE 
21  CONTINUE 

*******************************1c***4t-k*1c***1t**1t***1c**ic***1c**1ck**1c 
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11  CONTINUE 

ENDIF 
5000  END 

**************************************************** 
*  FUNCTION  SUBPROGRAM 

**************************************************************** 
DOUBLE  PRECISION  FUNCTION  DFUN2 (A, B, N, M, P, T, S) 
DOUBLE  PRECISION  A, B,  YT1 , YT2 
DOUBLE  PRECISION  P (90) , T (90) , S (90) 
IF     ( (A.LE.P(l) )    .AND.    (B.LE.T(l)))   GOTO  13 
IF     ((A.GE.P(N))    .AND.    (B.GE.T(N) ) )   GOTO  14 
IF     ( (A.LE.P(l) )    .AND.    (B.GE.T(M) ) )   GOTO  15 
IF     ((A.GE.P(N))    .AND.    (B.LE.T(N-M+1) ) )   GOTO  16 
IF     (A.LE.P(l))   GOTO  19 
IF     (A.GE.P(N))   GOTO  23 
MP  =  M+l 
GOTO  17 
23         I  =  (N-M+l) 
12=  N 
GOTO  22 

19  1  =  1 
12=  N 

22         DO  20  J=I,I2 

20  IF   (B.LE.T(J))   GOTO  21 

21  J  =J-1 

DFUN2=  S(J)+(B-T(J))/(T(J+1)-T(J) ) * (S ( J+l) -S (J) ) 
RETURN 

13  DFUN2  =  S(l) 
RETURN 

14  DFUN2=  S(N) 
RETURN 

15  DFUN2=S (M) 
RETURN 

16  DFUN2=S (N-M+l) 
RETURN 

17  DO  3  1=  MP,N,M 

IF   (A.LT.P(I) )   GOTO  4 

3  CONTINUE 

4  IF   (B.LT.T(I-M) )   GOTO  9 
IF   (B.GT.T(I-l) )   GOTO  11 
GOTO  12 

9         YT1=  S(I-M) 
YT2=  S(I) 
GOTO  18 

11  YT1=  S(I-l) 
YT2=  S(I+M-1) 
GOTO  18 

12  J1=I-M 
J2=I-1 

DO  5  J=J1,J2 

IF   (B.LT.T(J))   GOTO  6 

5  CONTINUE 


6  YT1=  S(J-l)   +(B-T(J-1) )/(T(J)-T(J-l) ) * (S ( J) -S ( J-l) ) 
I2=I+M-1 

DO  7  J=I, 12 

IF   (B.LT.T(J) )   GOTO  8 

7  CONTINUE 

8  YT2=  S(J-1)+(B-T(J-1) )/(T(J)-T(J-l) ) * (S ( J) -S ( J-l) ) 
18         DFUN2=  YT1+(A-P(I-1) )/ (P(I) -P(I-l) )*(YT2-YT1) 

RETURN 
END 

************************************************** 


*  SUBROUTINE  FOR  MATRIX  INVERSION 


SUBROUTINE  INVDET (C,N, DTNRM , DETM ) 

DOUBLE  PRECISION  C ( 11 , 11) , Jl ( 100) , PD, DD, CC, S , DETM, DTNRM 

PRINT*, '**C(1,2) ' ,C(1,2) 

PD=1.0D+00 

DO  124  L=1,N 

DD=0. 0D+00 

DO  123  K=1,N 

123  DD=DD+C(L,K) *C(L,K) 
DD=DSQRT (DD) 

124  PD  =  PD*DD 
DETM=1.0D+00 
DO  125  L=1,N 

125  J1(L+20)=L 
DO  144  L=1,N 
CC=0.0D+00 
M=L 

DO  135  K=L,N 

IF   ( (DABS (CC) -DABS (C(L,K) ) ) .GE.O.OD+00)   GOTO  135 

126  M=K 
CC=C(L,K) 

135  CONTINUE 

127  IF(L.EQ.M)   GOTO  138 

128  K=Jl(M+20) 
Jl(M+20)=Jl(L+20) 
J1(L+20)=K 
DO  137  K=1,N 
S=C(K,L) 
C(K,L)=C(K,M) 

137  C(K,M)=S 

138  C(L,L)=1. OD+00 
DETM=DETM*CC 
DO  139  M=1,N 
C(L,M)=C(L,M)/CC 

139  CONTINUE 
DO  142  M=1,N 
IF   (L.EQ.M)   GOTO  142 

129  CC=C(M,L) 
IF   (CC.EQ. 0. OD+00)   GOTO  142 

130  C(M,L)=0. OD+00 
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DO  141  K=1,N 

141  C(M,K)=C(M,K)-CC*C(L,K) 

142  CONTINUE 
144  CONTINUE 

DO  143  L=1,N 

IF   (Jl(L+20) .EQ.L)   GOTO  143 

131  M=L 

132  M=M+1 

IF   (Jl(M+20) .EQ.L)   GOTO  133 
136       IF   (N.GT.M)   GOTO  132 

133  Jl(M+20)=Jl(L+20) 
DO  163  K=1,N 
CC=C(L,K) 
C(L,K)=C(M,K) 

163  C(M,K)=CC 
J1(L+20)=L 

143  CONTINUE 
DETM=DABS (DETM) 
DTNRM=DETM/PD 
RETURN 

END 

****************************************************** 
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